Short-range correlations in quark matter 



F. Fromel^B and S. Leupold 1 - 2 

1 Institut fur Theoretische Physik, Universitat Giessen, Germany 
2 Gesellschaft fur Schwerionenforschung, Darmstadt, Germany 

We investigate the role of short-range correlations in quark matter within the framework of the 
SU(2) NJL model. Employing a next-to- leading order expansion in 1/N C for the quark self energy we 
construct a fully self-consistent model that is based on the relations between spectral functions and 
self energies. In contrast to the usual quasiparticle approximations we take the collisional broadening 
of the quark spectral function consequently into account. Mesons are dynamically generated in the 
fashion of a random phase approximation, using full in-medium propagators in the quark loops. The 
results are self-consistently fed back into the quark self energy. Calculations have been performed 
for finite chemical potentials at zero temperature. The short-range correlations do not only generate 
[ — , finite widths in the spectral functions but also have influence on the chiral phase transition. 
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To our present knowledge, quantum chromodynamics (QCD) [ij is the theory of the strong interaction. QCD 
is, however, not well suited for investigations at moderate energies. Due to the running coupling, perturbative 
calculations are not possible in this regime. Lattice QCD [2| is technically challenging and - in particular at finite 
chemical potentials - still limited in its applications [|[ . Consequently, effective models have been developed to explore 
the low energy phenomenology of QCD. These models treat the quarks usually as quasiparticles. Even in systems with 
■ considerable quark densities, short-range correlations that lead to a collisional broadening of the spectral function are 
f^*) \ not considered. 

The existence of short-range correlations in nuclear matter is well established. They are observed in A{e, f[p)X 
and A(e, e'pp)X experiments (see, e.g., [1, Q) and have been the subject of many theoretical approaches, see [sill for 
a review. The results of these studies agree rather well concerning the short-range effects. Even in a self-consistent 
calculation with a pointlike interaction [8|, [9( - using only one free parameter - good agreement with more realistic 
calculations and the experimental data is reached. Thus, what matters is not the detailed modeling but the overall 
strength of the interaction and the collisional phasespace. 

The success of using pointlike interactions in nuclear matter had motivated our first approach to short-r ang e interac- 
tions in quark matter [l(| on the basis of the Nambu-Jona-Lasinio (NJL) model. The NJL model [^[11, [III [SUB] * s 
. Z-t ', an effective quark model that respects the relevant symmetries of QCD. It is frequently used to investigate phenomena 
related to chiral symmetry. In the last years, it has also been used to explore the effects of color superconductivity - 
see, e.g., [l6l.[r7| for recent reviews. In the present work, as a first step, we ignore color superconductivity and explore 
the role of short-range correlations in quark matter in the chirally broken and restored phase. Clearly, for a complete 
picture also the phenomenon of color superconductivity - relevant at high chemical potentials /i and low temperatures 
T - should be taken into account. 

At zero temperature, just above the chiral phase transition, quark densities a few times larger than normal nuclear 
matter density are reached. It is not unlikely that short-range correlations exist at such densities and have influence 
on the medium [l8| . A simple estimate indicates that the character of the phase transition at low T changes when a 
small width is added to the quark propagator: The solid line in Fig. [1] shows the result of a Hartree-Fock calculation 
including a finite current quark mass. The 'S' shape of the curve is characteristic for a first-order phase transition 
[l9| . Somewhere in the shaded region the effective quark mass m* jumps discontinuously. The dashed curve has been 
obtained after inserting a constant width into the quasiparticle propagator, m* drops rapidly but continuously now 
- collisional broadening has turned the phase transition into a crossover. Note that Fig. [1] is meant as an illustrative 
exam ple. A more serious calculation with a self-consistently determined width will be given below. 

In [lfj] we have used the NJL interaction in a self-consistent calculation where the collisional self energy was 
determined by the Born diagrams, i.e., on a 2-loop level. We have found rather weak short-range effects in this 



'Electronic address: Frank.Froemel@theo.physik.uni-giessen.de 



2 



300 



200 



100 




320 



340 360 
H [MeV] 



380 



400 



FIG. 1: The quark mass as function of the chemical potential at T = in the NJL model. The solid line shows the Hartree-Fock 
result, the dashed line the result of a calculation with a constant quark width. The shaded area denotes a first-order phase 
transition. See the text for details. 



app roach - more than one order of magnitude smaller than in nuclear matter. A similar approach can be found in 
[20, HH, where the self energy was not determined self-consistently and the calculations were restricted to /i = 0. The 
results are comparable to [10(. 

Loop-expansions in the NJL model are questionable since the coupling is large. Expansions in the inverse number 
of colors, 1/N C , should be preferred 0, I n this work we bring our approach from [l(| to a consistent level in 
l/N c . Including all quark self energy diagrams of the next-to-leading order 0(1/N C ) in the self-consistent calculation 
leads to the introduction of dynamically generated mesons. In contrast to the standard Hartree+Random Phase 
Approximation (RPA) scheme [l3| , the mesons are part of the self-consistent procedure and their properties are fed 
back into the quark self energy. We restrict our numerical calculations to flavor symmetric quark matter at T = 
here. Nonetheless, the formalism can also be used at finite T. An approach to flavor asymmetric systems would be 
technically more involved but is theoretically feasible. See |24| for an approach to asymmetric nuclear matter. As 
already mentioned, diquark condensates and color superconductivity are presently also not considered. 

\/N r extensions to the NJL model on the quasiparticle level have been discussed in [22|,[23|, [25[. It has been shown 
23l l25j | that the generation of massless pions (Goldstone modes) in a perturbative expansion relies on a careful 



choice of diagrams for the meson polarizations. Our self-consistent ansatz iterates diagrams up to arbitrary orders 
in 1/N C but does not generate all contributions of any given order. Thus, the RPA pions become massive in the 
chirally broken phase. An approach that is self-consistent and yields massless pions (in the chiral limit) has not 
been proposed so far for the NJL model. The result of our self-consistent calculation could be used as input for an 
additional formalism which yields realistic pion properties. However, in the approach proposed in [23| and similarly 
in [2^, [27|, the Goldstone pions are not fed back into the calculation of the quark properties. 

In the present work we go beyond the quasiparticle models of [HI, [25| . We do not only consider the poles of the 
RPA spectral functions but also the continuum of off-shell states. As we will see later, these states are much more 
important for the quark properties than the poles. An analysis of the different on-shell contributions to the quark 
width will show that some of them are sensitive to the pion mass while others are not. Concerning the chirally broken 
phase, it will turn out that our self-consistent calculation yields a conservative estimate of the short-range correlations 
in quark matter. 

This work is structured as follows. We introduce our self-consistent 0(1/N C ) approach in Section[Tl] In Section Unl 
we show how to calculate the quark self energies and RPA polarizations using full in-medium propagators. The chiral 
properties of the RPA pions are also discussed. In Section IIVI we use quasiparticle approximations to explore the 
structure of the quark width qualitatively and to investigate the influence of the RPA pion mass. The technical 
details and the results of our numerical calculations - where we do not use any quasiparticle approximations - are 
presented in Section [V] A summary of our results and an outlook to future improvements are given in Section IVII 
The Appendices contain some of the more technical details of our approach. 
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TABLE I: Mean field NJL parameter sets with and without (chiral limit) finite current quark mass mo, using a three-momentum 
cutoff A The coupling G is fitted for Hartree calculations. The coupling Ghf for the Hartree-Fock approximation is 

obtained from G by a rescaling factor 12/13. m* and (uu) are the effective mass and the quark condensate that are found 
using the given parameters. 
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II. THE FORMALISM 
A. The NJL interaction 

The NJL model is by design the simplest effective quark interaction that resembles all relevant symmetries of QCD 
In the leading-order mean field approach - the Hartree+RPA scheme - a finite constituent quark mass 
is generated and breaks chiral symmetry dynamically at low /i and T, cf. Fig. [1] At higher \x and/or T a chiral phase 
transition occurs. The (Hartree+)RPA pions can be identified with Goldstone modes. 

The standard version of the two flavor NJL Lagrangian is given by 

£njl = i> (}$ - m ) i/> + G [(V^) 2 + (^ 75 r^) 2 ] , (1) 

where G is the constant qq coupling strength and the 73 are the isospin Pauli matrices, Too denotes a small current 
quark mass that breaks chiral symmetry explicitly. Due to the constant coupling, the model cannot be renormalizcd. 
Several regularization schemes exist [13j . the simplest one is the (three-)momentum cutoff A. 

G, A and too are free parameters. Their values are fixed such that reasonable results are obtained for the quark 
condensate and the pion decay constant in vacuum - and for the pion mass in the case Too 7^ 0. We list some typical 
mean field parameter sets in Table|H The simplicity of this model leads also to shortcomings like the lack of asymptotic 
freedom and the absence of confinement. Note that the cutoff is sometimes interpreted as a crude implementation of 
asymptotic freedom. 

The Lagrangian |T]) provides scalar and pseudoscalar qq interaction channels. It is also possible to construct a 
bosonized version of the NJL Lagrangian [lj, HH , 

£njl =ip(i0- Too) ip-ip (oTcr + 7r r + tt+P- + ip 

- ^ (a 3 + To + 27r +7 r_) , (2) 

where the meson fields are defined as a — 

-2G#VV and 7r 0,± = -2GipTo !± ip, and 

r CT = f a = 1 , r = f = ij5T 3 , T± = f T = ij 5 r± . (3) 

Note that t± = (n ± ir2)/v2- I<x,o,± are introduced here for further reference. They will be needed for the calculation 
of self energy and polarization diagrams. For a discussion on the quantization of |(2J) we refer to [28ll29j|. To investigate 
the phenomenon of color superconductivity fl7| , a vertex that introduces an attractive interaction in the qq channel 
must be added to |T]). We do not consider such an extension here. 

B. 1/N C expansion in the NJL model and RPA mesons 

The QCD inspired NJL model is a strongly interacting theory. GA 2 - the relevant quantity to estimate the 
interaction strength - has a value of approximately 2, cf. Table |TJ Thus, a perturbative expansion in terms of the 
coupling is not feasible. It is possible, however, to perform an expansion in the inverse number of colors, 1/N C . A 
simple scheme (see, e.g., [22|, [Hj]) allows to determine the order in 1/N C of diagrams by counting vertices and loops: 
A factor G is assigned to each vertex and a factor N c to each closed fermion loop. In this symbolic notation one sets 
0(G) = 0(1/N C ). The (free) quark propagator is of order unity. The order of a diagram with n vertices and to loops 
is then given by G n N™ = 0(N™~ n ), i.e., the difference of loops and vertices. 
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FIG. 2: The Hartree (a) and Fock (b) self energies and the direct (c) and exchange (d) diagrams of the Born self energy, (e) 
and (f) are diagrams of higher order in the coupling. The pointlike qq vertices have been replaced by finite interaction lines to 
show the structure of the diagrams unambiguously. 



FIG. 3: The 0(1/N C ) diagrams of the structure shown in Fig. [2] are added up in a random phase approximation to generate 
an effective meson exchange (dashed line). The first diagram is the Fock diagram, solid dots denote the NJL coupling (— 2G). 
The second line shows the RPA polarization 11;. 

Fig.[5]shows some diagrams that are of first (a,b), second (c,d), and higher orders (e,f) in the coupling. The Hartree 
diagram (a) is the only contribution of order 0(1) in l/N c to the quark self energy. The Fock diagram (b) as well as 
the direct Born diagram (c) are of order 0(1/N C ). The exchange Born diagram (d) is even of order 0(1/N^). 

The diagrams (e,f) are - like (b,c) - of order 0(1/N C ). The two diagrams were constructed by adding pairs of 
loops and vertices to (c). By adding further loops in (f), diagrams of any higher order in the coupling while still of 
order 0(1/N C ) can be found. The (D(l/N c ) diagrams that are constructed by dressing the quark lines in the 0(1/N C ) 
diagrams of Fig. [2] with the Hartree diagram, or by dressing the Hartree loop with the 0(1/N C ) diagrams will later be 
generated automatically in our self-consistent approach. They are not discussed explicitly to avoid double counting. 

Note that the full propagator of a self-consistent calculation can be dressed with self energies of arbitrary order. 
It remains 0(1) in leading order but gains contributions of higher orders in l/N c . In a self-consistent approach this 
means that one can select Feynman diagrams by their leading order to find a consistent set of diagrams for any 
given order. However, every diagram will contribute in subleading orders in a possibly incomplete way. The Hartree 
approximation constitutes an exception. Since the Hartree self energy is of order 0(1/ N®) the order of a (dressed) 
Hartree propagator remains exactly 0(1). 

Usually, all 0(1/N C ) diagrams of the structure shown in Fig. [U(b,c) and (e,f) are summed up in an RPA. They can 
be interpreted as contributions to an effective meson exchange as shown in Fig. [3] This allows us to treat a whole 
class of self energy diagrams in a simple and consistent way. The meson propagator is given by a geometric series, 



+ ... 





A, (A) 



2G [1 + (-2G) n,(fc) + (-2G) n ; (fe)(-2G) II,(fc) + . . .] 
2G 



i + 2Gn,(fc) ' 



(4) 



where n; (/ = a, 0±) is the polarization of the RPA mesons, cf. Fig. [31 Like the Lagrangian ||2J), the propagators have 
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FIG. 4: The set of Dyson-Schwinger equations corresponding to the fully self-consistent 0(1/N C ) approach. The double lines 
denote dressed quark propagators, dots denote the NJL coupling (— 2G). The second line shows the RPA like construction of 
the mesons (dashed lines), cf. Q. 

no kinetic part of the form k 2 — m 2 . The dynamics are governed by the underlying qq states and thus hidden in the 
polarizations. 

The propagators have poles at 1 + 2GRc 11/ (k) = 0. The poles correspond to bound qq states and are identified with 
the actual RPA mesons. At energies above ko = \i + m* (2m* in the vacuum), the propagators pick up contributions 
from unbound qq states. The mesons may decay into quark-antiquark pairs above this threshold. A wide range of 
states with large decay widths, the so-called qq continuum 28], is found. 

C. The 0(1/ N c ) approach 

The standard Hartree+RPA approach to the NJL model is of leading order in 1/N C . It can be summarized by the 
two Dyson-Schwinger equations of Fig. [¥J excluding the meson exchange diagram in the first line. The quark self 
energy is determined solely by the time-local Hartree diagram. Thus, the quarks remain quasiparticles and the RPA 
mesons have no influence on the properties of the quarks. Such a feedback can only be generated by a diagram of 
higher order in 1/N C . Note that the first Dyson-Schwinger equation - that determines the quark mass - is usually 
called gap equation in the Hartree approximation. 

In our investigation of the short-range correlations, we will determine the quark properties in next-to-leading order 
0(1/N C ). Therefore, the meson exchange diagram of Fig.[3]is included in the first Dyson-Schwinger equation in Fig.|U 
This equation will self-consistently generate all 0(1/N C ) contributions to the self energy. In addition, it introduces a 
feedback of the RPA mesons 1 - generated by the Dyson-Schwinger equation in the second line - into the quark self 
energy. The calculation of the RPA polarizations is part of the self-consistent procedure here. We will not use a pole 
approximation for the mesons since the qq continuum plays an important role in our approach. 

An analytic solution of the self-consistency problem presented in Fig. H] is not possible. A solution can only be 
found in an iterative, numerical calculation. By using full in-medium propagators - including off-shell contributions 
- we clearly go beyond the standard Hartree+RPA approach. The quark self energy has no longer the simple form 
of an effective mass. Due to the nonlocal meson exchange, it becomes complex and four-momentum dependent. The 
imaginary part of the self energy is generated by short-range correlations. It will be identified with the width of the 
quark spectral function. 

To keep our model numerically simple, we will use a three-momentum cutoff A [i~3j to regularize the divergent self 
energy and polarization integrals, cf. (|31I35I39[) . Like in [l(| and [2(| Hl[, we use the cutoff to restrict the three- 
momentum of each quark to values below A. A meson cutoff is implicitly included since the quark loop integrals for 
the polarizations are only finite for meson three-momenta below 2A. Note that the cutoff can also be implemented in 
other ways, see e.g. (23l. |25|. [26j . 

As long as the width does not become too large, most of the strength of the quark propagator remains in the vicinity 
of the on-shell peaks. The three-momentum cutoff serves then also as an energy cutoff, |po| < -^A < VA 2 + to* 2 . At 
higher energies, the on-shell peak will be located at three-momenta above the cutoff A. This is not a sharp cutoff 
when the peaks have a finite width, but the quark propagators at higher energies will be strongly suppressed in self 
energy and polarization integrals. 



On the Hartree level (for the quark propagators) the summation of all quark loop diagrams corresponds to a (classical) random phase 
approximation. We follow the literature, see e.g. [23l . and stay with the name RPA when propagators beyond the Hartree approximation 
are used — even when such RPA mesons do not satisfy the chiral theorems. 
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At this point, a remark concerning our earlier approach in [10] is in order. We have not considered dynamically 
generated mesons there. Instead, the direct 0(1/N C ) and the exchange 0(1/ N^) Born diagrams (cf. Fig. [2]) were used 
to construct the quar k self energy in a loop-expansion. In comparison to nuclear matter we found rather small widths 
for the quarks in [lCj. By including an infinite number of 0(1/N C ) diagrams in the new approach, the importance of 
collisional broadening will increase significantly. 



D. Quark fields and propagators 

We use the real-time formalism [3(| HH, HH, |33[ for our calculations. More details on in-medium techniques can be 
found in [34| and the transport theoretical approaches to quark matter [H, [5^, [35[ . In the following we assume that 
the system is in thermal equilibrium. For simplicity, we work only in the rest frame of the medium. 

The fundamental elements of our model are the causal, the anti-causal, and the non-ordered single-particle Green's 
functions, 

S^(l, 1') = -i (T[^(l)#(l')]) , 8^(1, 1') = -* (Va(l)^(l')) , 

S^(l, 1') = -i (T[^(l)^(l')]) , S<j(l, 1') = i • 

where a, (3 denote the spinor indices of the quark fields. Color and flavor indices have been suppressed here. T and T 
are the time-ordering and the anti-time-ordering operator, respectively. We also introduce the retarded and advanced 
propagators 

S ret (l, 1') = Q(h - iv)[S>(l, 1') - S<(1, 1')] = S c (l, 1') - S<(1, 1') , 

S av (l, 1') = Q(t v - t 1 )[S < (l, 1') - S > (1, 1')] = S c (l, 1') - S>(1, 1') . U 

Note that S rct and S av are complex conjugates of each other. 

For our calculations it will be advantageous to work in momentum space. We consider infinite systems in thermal 
equilibrium here. Hence, the Fourier transformed propagators will depend only on the energy and the momentum of 
the quarks. 

The relativistic quark propagators have a matrix structure in spinor space. The most general form of this Lorentz 
structure is found by a decomposition in terms of the 16 linear independent products of the Dirac 7 matrices (Clifford 
algebra). Considering the invariance under parity, time-reversal, and rotations, the Lorentz structure in the rest frame 
of the medium reduces to three independent components [36| , 

S(fc) = S s (k)t 4x4 + S (k) 7o ~ S v (k)e k ■ 7 = S s (k) + S^(k)^ , (7) 

with e*fc = k/\k \ and S l (k) = S v (k)e\.. S s , So, and S„ are scalar functions of fco and |fe|. In flavor symmetric quark 
matter, the functions do not have a structure in flavor space. Note that the sign convention in differs from [Toj | 
where S v had the opposite sign. S s , So, and S u can always be extracted from the full propagators S a /3 using 

trS(fc) tr[7 S(fc)] trjgfc ■ jS(k)} 
S«(k) = — 1 — , S (fc) = , S v (k) = . (8) 

Throughout this work we use the following convention: "tr" denotes a trace in spinor space while "Tr" denotes a full 
trace in color, flavor, and spinor space. 

The spectral function A(k) corresponds to the density of states at a given energy fco and momentum k. As the 
definitions in ([5]) show, iTr7 S > (l,l) and — iTr7 S < (1, 1) are related to the density of the free and the populated 
states, respectively [3l| . Thus, we can introduce the spectral function in terms of these propagators, 

A(k) = 1 [S> (k) - S< (k)] = -21m S rct (k) . (9) 

We use the definition ReF = (1/2) (F + 70^70) and ImF = (l/2i)(F - 70^70) her e [TSUI]. This definition ensures 
that ImS = (ImS s )l 4X 4 + (ImS^h^ 1 (i.e., 7/i = Re lfi ). 

The spectral function inherits the Lorentz structure of the Green's functions ([7]). The "normalization" of the purely 
real Lorentz components A Sl Aq and A^ can be found by using the methods in [371 ], 

—A (p) = l, / —A s (p)=0, / —A v (p) = Q. (10) 
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(13) 



Both poles of a retarded propagator are located in the same complex half plane [33j. When the integrand vanishes 
for large \po\, the integral can be closed along a contour in the other half plane. It can also be shown that 

f° tr S ret (p) = / f£ tr f , S ret (p)L, S r > - fc) = . (11) 
-00 27r J-00 27r L J 

The fermionic phasespace distribution function np(fc) determines the probability that a state is populated. In 
thermodynamical equilibrium, np(fc) depends only on ko, 

n F( fc o)- e(fco J /T + 1 — eo*-*D). (12) 

At T = 0, all states below are occupied while all states above fj, are free. The distribution function can be used to 
express S^(fe) in terms of the spectral function (30j . 

S>(k) = -iA(k)[l-ii F (k )}, 

S < (k) = iA(k)n F (ko) . 

In thermodynamical equilibrium, the retarded and advanced propagators in momentum space are given by 

^ret,av _ ^ 

orct.av/, \ _ 1 _ g + "^ rC ,aV (fc) /, .» 

M " # - mo - S rct ' av (fc) " fc r 2 cMv - m r 2 ct , av (fc) ' 1 ' 

where we have introduced the effective masses and four-momenta, 

m rct ' av (/fc) = m + Z T s et ' av (k) , fc^ ot ' av = k„ - £)f< av (/c) . (15) 

S rot ' av (fc) = E^ ct ' av (fc) + Hjf ' av (fc)7 M are the (complex) retarded and the advanced self energy, respectively. They 
have the same Lorentz structure as the propagators ([7]). In contrast to the time-ordered propagators S c,a that gain 
extra terms (33[, S ret,av keep the simple form (|14|) in the presence of a medium. 

For our calculations it will be necessary to disentangle the real and imaginary parts of the propagators, e.g., to 
determine the spectral function. The real and the imaginary part of the denominator on the rhs. of Eq. (TTi|) can be 
identified with 

V{k) = ReP - Rem(fc) 2 - [lmS^ et (fc)ImE^ t (A;) - ImSf(fc) 2 ] , (16) 

WO) = -2 [Refc M Im ££,<.(>) + Rem(k)ImY, r s et (k)] (17) 

for the retarded propagator. Using Y,* ct (k) = E av (£;), we find V av (k) = V{k), W av (k) = ~W{k) for the advanced 
propagator. Note that Rem = ReTO rct ' av and Kek = Refc ret,av . The expansion with the complex conjugates of the 
denominators yields [H| 



(% 1Ct aV + m rot - av (fc)) (V(k) =F i W(k)) 



S ret > av (fc) = ^ = ' \ — . (18) 

The upper sign in the numerator refers to the retarded and the lower sign to the advanced propagator, respectively. 
In both cases the functions V and W of ()16|17|) are used. 

The zeros of V(k) determine the poles of the propagators. W(fc) generates the width. In the following, we will refer 
to the broadened poles as on-shell peaks. The on-shell energy k^ s (k) and the on-shell momentum k os (k ) are defined 
by 

=0, 7>(fc ,£os(fco)) =0- (19) 

We note that V and W both depend on the full self energy, not only the real or imaginary part. As we will see 
later, Im E rot does not become too large in the on-shell regions. Hence, the influence of the quadratic terms in (fl6|) 
is limited. 

The real and the imaginary part of the propagators can be easily extracted from Eq. (|18p . We identify the Lorentz 
components of the spectral function with 

v ' -p 2 + yy 2 ^ y ' p 2 +VV 2 y ' 

In the absence of a current quark mass, the Lorentz scalar components S s and E s vanish when chiral symmetry is 
restored. The Lorentz structure is then simplified even further. In case of a finite mo, S s and E s will become small 
but do not vanish entirely. 
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E. Meson propagators and spectral functions 

The a and tt 0,± Green's functions in terms of the meson fields [28| can be defined in analogy to the quark Green's 
functions (|5l6p . Like for the quarks the a and spectral functions are introduced via the non-ordered 

propagators A^, 

- i [A>(fc) - A<(fc)] = -21m A rct (fc) . (21) 

We have dropped the indices a, 0, ± since the definition is the same for all mesons. Note that A rct * = A av . In 
thermodynamical equilibrium, the bosonic distribution function 

n B(fco)= eW r_i ^ -e(-fco), (22) 

is - like np - a simple step function at zero temperature. We can use n B to express the non-ordered propagators A^ 
in terms of the spectral functions [3(| , 



A>(A) = -^(fc)[l+n B (Ao)] J 
A<(fc) = -i^(fe)n. B (fto) . 

Note the subtle differences to the relations for the quarks. The distributions np in (fT3| correspond to — n B in (|23p . 
i.e., Pauli blocking is replaced by Bose enhancement. 

In Section llI B| the RPA meson propagators have been determined in terms of the NJL coupling and the polarizations 
n, cf. Eq. Q, 

9C 

AT"* (A) = — t , (24) 

1 y ' i + 2Gn[ ot (fc) v ' 

where I = a, 0, ± denotes the type of the meson. We can easily disentangle the real and the imaginary part of the 
retarded propagators. The spectral function £ = — 2ImA Ict is identified with 

«»)_ . (25) 

(jij + Rcnr'M) +imnr'M 2 

The term ^ + ReIIJ ct determines the pole structure of £/. The imaginary parts of the polarizations generate the 
width. The meson spectral functions are not normalized to a certain value. There are no explicit coupling factors 
that enter the quark-meson interaction vertices, cf. Eq. @. Hence, the RPA spectral functions can be interpreted as 
the products of normalized spectral functions and the effective quark-meson couplings. 



F. Self energy, polarizations and widths 

In analogy to the single-particle Green's functions (0, we can define causal, anti-causal, and collisional (non- 
ordered) self energies [3l|, HH, [33[ . We have already encountered the retarded and advanced self energies £ rct ' av . 
They are defined as 



S ™t(l ) l') = S nrf (1}1 / )+e(tl _ tl#) |- S > (1)1 / ) _ E < (1)1 / ) ] 

E av (l, 1') = S mf (l, 1') + G(t v - h) [£<(!, 1') - £> (1, 10] 



In our 0(l/N c ) approach, the mean field self energy S mf is identified with the Hartree self energy E H , cf. Fig. [U It is 
mainly responsible for the generation of a constituent quark mass. The collisional self energies correspond to the 
meson exchange diagram. They generate the short-range effects that arise from decays and collisions in the medium. 

The Fourier transformed E ret (A:) enters the spectral function A(k) = -2ImS rct (fc) pJ]). As wc can see in (|17I18|I . 
the width of the spectral function is determined by — 2ImE ret (/c). In momentum space we find — 2ImE let (/s) = 
i E > (/c) — i E < (fc) [31(, a relation similar to ([9]). Therefore, we define the quark width T by the sum of the collisional 
self energies, 

r(jfe) =iH > (k) -iE K (k) = -2ImE rot (fc) . (27) 
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Note that — iS < (fc) and iE > (fc) are identical to the total collision rates for scattering into (gain rate) and out of 
(loss rate) the configuration (k ,k), respectively [l(| H(j. The relation on the rhs. of (|2"T)) is just the optical theorem 
[l|, [HI]. It relates the imaginary part of a Feynman amplitude (ImE rct ) to a total cross section or collision rate 
(«S > -iS < ). Cutting the meson exchange diagram in Fig. [3] yields processes like q ^ qn and qq — ► cr. We will 
examine these processes in more detail in Section IIVI 

Like the quark propagators ([7]) and self energies, the width consists of three Lorentz components T Si o tV which are 
functions of fc and \k |. Using the definition of the on-shell energy and momentum (|19p . we can introduce on-shell 
self energies and an on-shell width, 

Sos(fco) = E(fc 0) ^os(fco)) , r os (fc ) = r(fe , fcos(fco)) , - 2g . 
s os (fc) = s(fcg s (fc),fc), r os (fc) = r(A : ° s (fc),fc). 

The RPA polarizations do not include a time-local contribution. Thus, the retarded polarization IIJ ct can be 
expressed solely in terms of the collisional polarizations IL^ , 

n rot (i,T) = e(< 1 -< 1 0[n>(i,i')-n<(i,i')] . (29) 

Like S^, the collisional polarizations zll > and —ill < can be identified with total collision rates. We find processes 
like 7r — ► qq, aq — > q, etc., when cutting the polarization diagram in Fig.O Thus, we introduce the RPA meson widths 
by 

r,(*) = [iu>(k) - in<(fc)]/(2fe ) = -imn[ ct (fc)/fc . (30) 

The additional factor 1/fco in comparison to (|27[) is required to get the dimension of the widths right. The factor 1/2 
is conventional. It ensures that the nonrelativistic limit of this definition is identical to the widths from nonrelativistic 
models III EL 



III. CALCULATION OF SELF ENERGIES AND POLARIZATIONS 



A. Mean field self energy and quark width 

In our 0(1/N C ) approach of Fig. [4] the Hartree diagram takes the role of the mean field self energy E mf . It is 
calculated using a full in-medium propagator that includes a finite width, 

-iE mf =-iE H = 2G J ^TTS c {p) = 8iGN f N c J A s ( P )n F (po) ■ (31) 

Due to the trace, only the isospin scalar (I = er) NJL vertex contributes to E H . On the rhs. we have used Eqs. ()6I7|11|) . 
The Lorentz structure of E is simple - it has only a Lorentz scalar component. Note that the integral must be 
regularized by a cutoff. 

The definitions of the quark condensate and the quark density are closely related to E mf . For the quark condensate 
we find in isospin symmetric quark matter, cf. (|31|) . 

(««) = (dd) = \ (H>) = / TrSC (p) = -4G" sH • (32) 

The quark number density p([i,T) is related to (V'^V'a) = —iT^7o S < (1, 1). Using (fLl)) the density can be calculated 
by integrating the spectral function weighted with the occupation propability over all states with a positive effective 
energy Repo $15]), 

p(M) = 4A/iV c / /°° ^A (p)n F (p ) . (33) 

J A JRep a =0 27T 

Aq is the only Lorentz component of A that remains finite in the nonrelativistic limit. In a quasiparticle ap- 
proximation, the peaks of Aq become ^-functions, see (|45|) . At zero temperature the integral can then be solved 
analytically, 

Pqp (fi) = N f N c /(3n 2 )k 3 F , (34) 

where kp = \/ p 2 — m* 2 is the quasiparticle Fermi momentum. Note that the Fermi momentum of the 0(1/N C ) 
approach is defined by V(jjl, kp) = 0, cf. (|16I19|) . 
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B. Collisional self energies and quark widths 

The collisional self energies are calculated from the meson exchange diagram in Fig. [3l using full in-medium 
propagators. The Feynman rules of the real-time formalism [3l| yield 

-iE^(fc) = £ J ( ^ 4 {2v)H\k -p + r)r ; S^ (p)f, Af (r) . (35) 

The contributions from the different RPA pions are equivalent up to the different propagators Aj^, A^ and the isospin 
factors in T;, T; In flavor symmetric quark matter, the isospin matrices t; commute with the quark propagators. 
Moreover, we have Ao = A + = A_ since the polarizations IIo,± are equal, cf. I|40p . Thus, it is not necessary to 
treat the different pions separately here. After introducing universal pion propagators A<, their contributions can be 
added up, 

-iE*(A) = J ^[sf(p){Af(p-k)-3Af( P ~k)} 

+ S$( P )r{Af(p-k)+3Af(p~k)}] . (36) 

To determine Im S rct and the quark width, it is convenient to replace the propagators in (|36|) by spectral functions 
according to Eqs. (|13I23|) . The width is given by the sum of the collisional self energies as shown in (j27|) , 

d 4 p 



r «./*(*) = I T^ZTjMPo) + Mpo ~ ko)}A s Ap) [&(p - k) T 3^(p - k)] , (37) 



where the upper sign refers to the Lorentz scalar and the lower sign to the Lorentz vector component. Ai{p) should 
be read as A v (p) cosde], when determining IV 

For T = 0, np and ne become simple step functions (| 12I22[) . This can be used to combine the distribution functions 
with the limits of the p integration, 



dpo[n F (p ) + n B (p - k )]- ■ ■ / dp - ■ ■ ■ (38) 



ko 



The width is zero at ko = fJ. in this case: States at the Fermi energy are stable (quasiparticles) . Due to Pauli blocking 
it is not possible to scatter into our out of those states. For finite temperature this restriction is less strict. At low 
temperatures, however, the width will remain small for energies in the vicinity of the chemical potential. 



C. Collisional polarizations and RPA meson widths 



For the collisional polarizations Hf(k), cf. Fig. [3) we find 



■iU<(k) = - 



d A p d 4 q 
(2^)4 (2^)4 



(27r) 4 (P(fc - p + q) Tr %T t i (p)iTii S>(g) 



(39) 



Working out the trace yields in flavor symmetric quark matter, where the isospin matrices r; commute with the quark 
propagators, 



-iUfjk) = -4N f N e 



d 4 p 



± Sf (p) Sf( P - k) + S<(p) S»(p - k) 



(40) 



The upper sign in the integrand refers to a and the lower sign to the 7r case. Since t^t± — 1 ± T3 and Trr3 = 0, the 
same result is found for all pions. 

From Eq. (H|) and IT^ = IT^E = II f follows that Ag-(fc) = A^(k) = A<(fc). Thus, we can introduce a universal pion 
propagator A n , as already done in (|36p . This simplification is not possible in flavor asymmetric matter. There, the 
quark propagators have a nontrivial structure in flavor space and will not commute with t\ . 



11 



The widths of the RPA mesons are found with the help of Eq. (|30j) . We replace the quark propagators in (|40f by 
spectral functions and distribution functions n F , 

„. 2N f N c f d 4 p , . , . , 
rV(fc) = ^ y ^4 [n F (po -ko)- n F {p )\ 

x [±A(p)A(p-fc)+^(p)^b-fc)] ■ (41) 

Note that the integral is antisymmetric in k. To see that, one has to perform the substitution p — * p— ^k. n F becomes 
a simple step function (|12[) at zero temperature. This can be used to rewrite the po integration, 

+oo t ^ r^+h, 

dpo[n F (p - k ) - n F (p )]- • • * / dp - ■ ■ . (42) 



D. Dispersion integrals 



In the NJL model, cutoffs are used to regularize the self energy and polarization integrals. Such cutoffs yield a 
violation of analyticity when we calculate both, the real and the imaginary parts of S rct and n[ et , using Feynman 
rules. Some properties of the spectral functions - like the normalization of A(k) (fTOjl - would be lost. 

Our main interest in the present work is to investigate the collisional broadening of the spectral functions. For 
the self-consistent approach it is more reasonable to preserve analyticity than to calculate the real parts directly. 
Therefore, we will calculate only Im£ lct (r SjA1 ) and Imil) ct from the Feynman rules as shown in Sections IIII Bl 

and llll Cl The real parts will be determined by the dispersion relations 

ReS rct (fc ,fc) = S mf +— V [ +0 ° dp T } P0,k ^ + const. , 
2tt k -p 

o met/, T\ 1 t> [°° , 2Po r iiPoik) 

ReU, (k ,k) = -V / dpr* — -5 5 h const. 

Jo k o- Po 

In the second relation, we have used the antisymmetry of poTi that has been mentioned after Eq. (|4ip . Constant 
- or at least energy independent - contributions to the real parts do not interfere with analyticity and cannot be 
calculated from dispersion integrals. Thus, the (purely real) mean field self energy must be calculated separately. 

As indicated in (|4"5)) . there exist further energy independent contributions to ReS ret and ReII[ ct . They should not 
be ignored. For example, the RPA pion mass depends strongly on the inclusion of such shifts. Identifying constant 
terms - like the Fock diagram that is part of the meson exchange in Fig. [3]- requires an investigation of the real parts 
that are given by the Feynman rules. The details of this rather technical exercise can be found in the Appendices [B] 
and[C] The full dispersion relations including all energy independent shifts are given in ([C6]) for Re E ret and in (jB22)) 
for Renf et . 



E. Masses of the RPA mesons 



As discussed in Section III Bl we use the RPA mesons mainly as a tool to generate the next-to-leading order quark- 
quark interactions. Nonetheless, it is interesting to explore the properties of the mesons. The quasiparticle approaches 
of [13, [2^,[2^| have shown that the RPA pions, i.e. the bound qq states, may loose their Goldstone boson character in 
1/N C extensions of the NJL model. In Appendix [A] that is briefly summarized below, we demonstrate how this effect 
arises on the mean field level. 

The poles of the retarded RPA propagators are given by the zeros of l + 2GRenf ct (fc), cf. (21|). In the Hartree+RPA 
approximation (see (|A4IA7p in Appendix |A1 below) . the term 1 + 2GII n = (m* — S mf )/m* becomes zero in the chiral 
limit (m = 0), cf. (KE\i . The poles are determined by Re IL™^ ~ fc 2 lml(fc) and Re 11^ ~ (fc 2 - 4m* 2 )ImI(A:). Hence, 
the RPA pions are massless and the RPA sigma has a mass of 2m*. For finite mo, 1 + 2GII n will generate only a 
small shift and the RPA pions remain light. The effect is shown in Fig.O 

In the 0(1/N C ) approach, ReII[ ot is given by (|B22[) and (|BlllB2ip . As numerical checks confirm, the properties of 
H-^tt ~ [fc 2 I 10t — Jct°tt] do not change too much. J™* remains very small and J™* is approximately given by 4M 2 I ret (fc), 
where M = rn + ReS™' s (fc = 0) is the constant (real) part of M from (|B14p . However, the term 1 + 2Gn n will - 
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FIG. 5: The real part of the denominator of the RPA pion propagator in the Hartree+RPA and Hartree-Fock+RPA approxi- 
mations, using parameter set I from Table [I] The zeros determine m^, see the text (and Appendix [XJ for details. 



like in the Hartree-Fock+RPA approach (|A6|) - not vanish for mo = 0. Eq. (|B11[) yields 

™o + S const . - £ 



l + 2GIT n 



m + ^const. 

16GNfN c f d 4 p £^ ct (p)-£ 



Im a - 2 const n F (po). (44) 



m o + Sc 0ns t. J (2tt) 4 ' p 2 ct - m 

The first term corresponds to the mean field result (|A4IA7[) with m* — m + S^ onst . The second term appears only 
in the full calculation, where E ret becomes complex and four-momentum dependent. 

It is reasonable to chose Sconst as large as possible (see the discussion below (jB10p ) to minimize the second term. 
Due to next-to-leading order contributions, S^ onst will then not be equal to E m . Using a numerical result of our 
approach, S mf w 0.75 S^ onst , leads to values of 0.15 — 0.2 for the first term while the second term remains small. 
Hence, the RPA pions are massive even in the chiral limit. Qualitatively, this resembles the Hartree-Fock+RPA 
approximation discussed in Appendix [A] see Fig. Quantitatively, the shift exceeds the Hartree-Fock+RPA value 
of 0.08. 

The large shift raises the question which influence the RPA pion mass has on the quark properties in our self- 
consistent calculation. We will address this problem in Section ITVl There, it will be found that the off-shell states 
of the pion propagator - and not the bound qq states - have the largest influence on the quark self energy. The qq 
continuum states in the timelike region are not generated by poles of the pion propagator but by large decay widths. 
They are largely independent of the RPA pion mass and differ within reasonable limits in the 0(1/N C ) and the 
Hartree+RPA approaches. The large RPA pion mass may lead, however, to a suppression of the contributions from 
the spacelike region of the propagator to the quark width. Our analysis in Section IIVI indicates that our approach 
will - on average - underestimate the short-range effects in the chirally broken phase. 

Let us conclude this section b y re viewing why the chiral properties of our model are broken and how they may be 
restored. It has been shown in [23L [25[ that the chiral properties of dynamically generated pions in a perturbative 
expansion depend on a careful choice of polarization diagrams. Only when certain contributions cancel each other, 
the pions become massless in the chiral limit. In 1/N C extended schemes, the quark propagators are not purely of 
order 0(1), cf. Section I'll B I When the RPA polarizations are calculated using these propagators, contributions in 
subleading orders will be automatically generated - in an incomplete way. Thus, the cancellation effect also remains 
incomplete beyond leading order and the RPA pions become massive. 

A 1/N C extension to the NJL model that restores the chiral theorems in next-lo- leading order can be found in 23]. 
The (quasiparticle) approach is based on Dyson-Schwinger equations similar to those in Fig. H] (see Figs. 1 and 2 
in [23|). These equations, that yield massive RPA pions, are solved self-consistently as a first step. The physical 
mesons are found in a second step, where the 1/N C extended polarizations shown in Fig. [B]- including next-to- leading 
order corrections 2 - are calculated using the self-consistent quark and RPA propagators. The 1/N C improved pions 



2 Note that a restoration of the chiral properties in higher orders may be possible by using Bethe-Salpeter equations that generate vertex 
corrections - like those in Fig.(6fc,d) — in all orders |27| . 
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FIG. 6: The 1/N C extended polarizations of [231 ] . Double solid and single dashed lines indicate self-consistently calculated quark 
and RPA meson propagators, respectively. Double dashed lines indicate 1/N C extended mesons. See the text for details. 

are not further modified, i.e., they are not reinserted into a self-consistent calculation. This would destroy the chiral 
properties again. 

Following the approach of [23[ , we could determine the Goldstone modes of our model (the diagram in Fig. [U(b) is 
generated automatically). However, this would have no influence on the properties of the quarks - the main interest 
of the present work. They would still be calculated from the Dyson-Schwinger equations of Fig. [5] where the massive 
RPA pions enter. Hence, such an extension is of limited use for us. 

A fully self-consistent approach that yields massless pions has not yet been proposed for the NJL model. (In [4lJ 
such a scheme has been proposed which, however, does not work in the vicinity of the phase transition.) Adding the 
vertex corrections of Fig. [S](c,d) to the second Dyson-Schwinger equation in Fig. 0] should improve - but not fix - the 
chiral properties of the pions and generate a feedback on the quarks. However, the numerical efforts would increase 
significantly by involving two- and three-loop diagrams that must be calculated in every iteration. Since the numerical 
calculations are already rather involved we do not pursue this strategy any further at the present stage. Accepting 
that our approach yields a conservative estimate of the short-range correlations in the chirally broken phase, we use 
the (uncorrected) RPA mesons in our calculations and ignore the vertex corrections for the time being. 



The self energies — i S < (fc) and i S > (A;) are identical to the total collision rates for scattering into and out of a quark 
state with four- momentum k [3(3]. Due to the variety of contributing scattering and decay processes, - and thus 
the width F = i S > (fc) — i E < (k) - have a complicated structure and a non-trivial energy dependence [42[ . For a better 
understanding of this structure, we will investigate the collision rates on the level of the Hartree+RPA approximation 
at zero temperature. The results will be helpful later when the results of the full calculation are discussed. In addition, 
they allows us to estimate the influence of a too large RPA pion mass on the quark properties. 



peak at ko = E^ and an antiquark peak at k = Ef.. An antiquark corresponds to a hole in the populated quark 
states. Since the chemical potential is positive in our calculations, there are no holes in the Dirac sea. Consequently, 
there are no antiquarks in the medium at T = 0. 

Each of the meson spectral functions ^ a ^ n f|25|) consists of two components that we will treat separately. Bound qq 
states generate the poles of the RPA propagators and show up as peaks in £ CT . X . In addition, the spectral functions 
include broad ranges of unbound qq states, cf. Section [II Bt The qq continuum is located at |fco| > fi + m*, i.e., above 
the threshold for the decay into qq pairs (28|. Landau damping (nq — > q, etc.) generates further contributions at 



IV. QUARK AND MESON SCATTERING 



A. Mean field spectral functions 



In the Hartree approximation, the quark spectral function A consists of two quasiparticle peaks, 
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spacelike four-momenta [43} . In those regions, the spectral functions are characterized by large widths r CT 7r and not 
by poles. 

The RPA pion peaks are located below the qq decay threshold. Hence, they generate pronounced contributions to 
the collision rates that are well separated from the continuum contributions. For the quasiparticle component of the 
RPA pion spectral function we can use an expression similar to (I45p . 

Uk)~6(k -E£)-6(k +EZ), (46) 

where E k = (fc 2 + m 2 ) 1//2 . is the effective pion mass that is determined by the poles of the RPA propagator. We 
will refer to a meson with negative energy fco as a meson with positive energy fco = — fco in the following. Since we 
are only interested in energy thresholds, it is not necessary to specify £ w in more detail. To investigate the off-shell 
contributions to the quark width, we will later replace the incoming and outgoing mesons in (|35() by qq pairs. The RPA 
sigma peaks are located above the qq decay threshold and thus are rather broad. For our qualitative investigation, it 
is sufficient to investigate the qq continuum contributions in this case. 

B. Thresholds 

The collisional self energies in Eq. ([55)) can be immediately identified with total collision rates: The propagators 
—i S < and i S > (fT3"[) are just the densities of the populated and the free quark states. The meson propagators can be 
interpreted accordingly, cf. (|23[) . Hence, the loss rate iT, > (k) corresponds to a quark with four-momentum fc that is 
added to the system and scatters off a meson into a free quark state. The gain rate — zS < (fc), on the other hand, 
corresponds to decays of a quark from the medium into a quark with four-momentum k and a meson. 

Note that incoming quarks with negative energy correspond to outgoing antiquarks (and vice versa) [421 ] . The 
(5-function in ([35]) yields fco + ro = pa, corresponding to the processes kr p. Since the energy of an antiquark is 
given by ko = —fco, the condition should be reinterpreted as ro = fco + po when fco is negative. This corresponds to 
the processes r kp. Accordingly, an incoming meson with negative ro is reinterpreted as an outgoing meson with 
positive energy ro = — ro. 

To investigate the contributions from the RPA pion peaks to the collision rates, we replace the propagators in ([3"5[) 
by the spectral functions ([45I46[) and the appropriate distribution functions, cf. ([13l23p . This generates four terms 
with different combinations of the quasiparticle peaks for each of the two self energies S^. At zero temperature, some 
terms are ruled out by the combination of <5- and step functions in the propagators. In total we find 

E p - E;_ k ) [1 - n F (E p )} n B (E p - fc ) , (47) 

E p + E;_ k )n F {E p ) [1 + n B (E p - fc )] 

+ / d 3 p S(k + E p + E;_ k ) [1 + n B (-E p - fc )] . (48) 

J A 

The corresponding processes are shown diagrammatically in Fig. [JJ There is one process (a) that contributes to the 
loss rate i S > and - taking into account that fco can be positive or negative in the first term of ([48[) - three processes 
(b-d) that contribute to the gain rate -iE < (antiquark loss rate). In the last term of (|4"8"[) . the Bose distribution 
demands fco < — E p . Hence, the integral will only be finite for negative fco- 

We can now determine the energy thresholds of the processes in Fig. [JJ Besides four-momentum conservation, 
the integrands of (|47I48[) include the following constraints: Incoming quarks from the medium must have a positive 
energy below the chemical potential. Outgoing quarks must have an energy higher than [i (Pauli blocking) . Outgoing 
antiquarks and mesons must have an energy that exceeds their mass. In addition, the (anti-)quark three-momenta are 
regularized by the cutoff A of the model. Since (j45[) provides a sharp relation between the energy and the momentum 
of an (anti-)quark, the energies po and po are restricted to values below E\ = V A 2 + m* 2 . 

In Table [TT] we summarize the thresholds that can be found for the processes of Fig. [Jj Where two values are given 
for the same threshold, the first one follows from Pauli blocking or the upper energy limit for quarks from the medium 
and the second one is a kinematical constraint. These limits are complementary. Which one is stricter depends on |fc |. 
The thresholds in Table HT1 are not always the strictest limits that can be found. However, to identify the contributions 
from individual processes to the quark width they are perfectly suitable. 

The continuum of off-shell states in the RPA spectral functions cannot be represented by quasiparticles. To 
investigate their contributions to the collision rates in a similar way as the contributions from the bound states, we 



zS>(fc) 



d 3 p <5(fc - 



-iZ<{k) ~ ( d 3 p S(k 



15 




FIG. 7: Scattering and decay processes corresponding to (k) (|35[) . The thick lines carry the external four-momentum k. See 
the text for details. 



TABLE II: Energy thresholds of the processes in Fig. see the text for details. Note that E\_ k = [(A + \k |) 2 + raj] 1 ' 2 , where 
m,r is the pion mass that is determined by the poles of the RPA propagator. 



Process 



Thresholds 



(a) 



k 2 + (m* + m-n 



< 



hi 



E A + El_ 



(b),(c) 
(m* > rriir) 



At - J(\k\+k F ) 2 +m 2 



h) 



fl — 



^Jk 2 + (m* - m^) 2 



(b),(c) 
(m* < m,r) 

(d) 



\k | - k F ) 2 +m 2 — n 



k 2 + (m* + niir 



E 



< 



A+k F — A* 



\Jk 2 + (m* — m-n) 2 
Ea + E\_ h 



rewrite the integrals in Eq. I|35p . Thereby we replace the outgoing mesons by (quasiparticle) qq pairs. Using the 
relation Af = A[ ct uf Af v HHEI in Eq. j35|) we find 



d i p d q d^r 
(^r) 1 (2^)4 (2^)4 



(27r) 4 (5 4 (fc+p- q-r) 



(49) 



xr ; S^(r)f ; A[ ct (r-fc)Tr f , S%)r z S*(«) Af v (r-fc). 



Diagrammatically, the self energies (|49|) have the form shown in Fig.[8j Eq. (|49|) and the diagrams in Fig. [8] resemble the 
approach of [ioj ]. where we have calculated the collisional self energies from the Born diagrams in Fig. [5] Here, however, 

the bare NJL coupling G 2 is replaced by the propagators of the dynamically generated mesons, A rct A av = | A rct f . 

The thresholds are now generated by the incoming and outgoing (anti-)quarks that are explicitly given by the 
propagators S^. Hence, we can ignore the off-shell structure of the meson propagators in (|49|) . It is also not necessary 
to distinguish between a and n contributions - they differ only in magnitude but not in the thresholds. 

Replacing the quark propagators in (|49[) by mean field spectral functions (|45[) and Fermi distributions, we find the 
processes shown in Fig. [9l Many of the formally possible processes are ruled out by the distribution functions (no 
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FIG. 8: Self energy diagrams corresponding to — iE > (left) and — i E < (right) in the rewritten form of Eq. (|49 
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FIG. 9: Scattering and decay processes corresponding to S^(fc) in the form of Eq. (|49|l . Thick lines carry the external 
four-momentum k. See the text for details. 

incoming antiquarks) . For simplicity, we do not display u-channcl processes in Fig. [5] They have the same thresholds 
as the shown t-channel processes (the full calculation includes both channels). 

The processes (a) and (b) in Fig.[S]are contributions to the loss rate i S > (fc). (c) and (d) are two of the contributions 
to the gain rate — iE < . Note that they differ only in the sign of feo- The processes (e)-(g) are the remaining 
contributions to — i S < for negative fco- The thresholds for these processes are found by applying the same constraints 
to the incoming and outgoing quarks as in the analysis of the processes in Fig. [Jj We have summarized all thresholds 
in Table [ml 



C. Density dependence 

Eqs. ([55)1 and allow us to make some estimates on the density dependence of the processes in Figs. [7] and [5] 
Each incoming quark from the medium can be identified with an integral over S < for positive energies, cf. (|33p . and 
thus introduces a linear density dependence into the corresponding collision rate. Each outgoing quark - not counting 
quarks with the external four-momentum k - corresponds to an integral over S > , i.e., the density of the free states. 
Due to Pauli blocking and the (fixed) cutoff, this density decreases when fi rises. 

It follows from these considerations that the growth of the collision rates - and the short-range correlations - will 
be limited. We have already investigated this effect for nucleons in nuclear matter [46|. There, Pauli blocking leads 
to a saturation of the correlations at a density of 2 — 3 times normal nuclear matter density. 

Pauli blocking becomes effective as soon as a medium is present. Its influence increases continuously as /j, rises. 
The cutoff acts only on outgoing states with large energy. When the total incoming energy is low, such states cannot 
be reached. Depending on the considered process, the cutoff may be only relevant if fj, and/or fco axe large. 
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TABLE III: Energy thresholds of the processes in Fig. [9j see the text for details. 



Process Thresholds 



(a) fi < ko < 2Ea — m* 

2/i + m* I 

(b) \ < k < 3E A 

y/k 2 + (3m*) 2 J 

(c) ,(d) 2/i- \f(\k\ + 2k F ) 2 + m* 2 < k < /J, 
(e),(f) m* k 2E A - m* 

(i + 2m* 

(g) . } k Q 3E A 



^k 2 + (3m*) 2 



We use process (c) of Fig. [9] as an example to illustrate the role of the cutoff: In this process the cutoff becomes 
relevant for p > E\, i.e., when the incoming energy exceeds £a + k . The maximum energy of the two incoming 
quarks is 2/i. Hence, the cutoff suppresses only contributions to the collision rates when 2/i > E\ + fco (> E\). For 
lower /z, final states with energies above E\ cannot be reached kinematically. 

Note that E\ = (A 2 +m* 2 ) 1 / 2 is - due to to* - density dependent and jumps at the phase transition. A numerical 
check shows that /i w E\/2 in the chirally broken phase. Above the phase transition, however, /i is well beyond this 
limit. Therefore, we can expect that the density dependence of process (c) in Fig. [HI and some of the other processes 
changes significantly at the phase transition. All these considerations will help us to interpret our numerical results 
presented below in Section IVl 



D. On-shell width 



The on-shell peaks dominate the structure of the quark spectral function. Hence, processes that contribute to the 
on-shell width have the largest influence on the properties of the medium. This is of particular interest for the present 
0(1/N C ) approach: The thresholds of the processes involving bound qq states depend on m^, cf. Table Hfl As we 
already know, the mass of the RPA pions in our approach will be larger than that of the physical pions in the chirally 
broken phase - and thus shift the thresholds. Furthermore, we will see below that the magnitude of the contributions 
from some of the processes in Fig. [S] depends on the pion mass. This will have some impact on the on-shell width and 
the properties of the medium. 

A process contributes to the on-shell width when the on-shell energy Ek = (k 2 + to* 2 ) 1 / 2 is located within its 
thresholds. As Table HTT1 shows, the processes (a), (c), and (e,f) of Fig. [5] contribute to the on-shell width of all free 
quark states, all populated quark states, and all antiquark states, respectively. The processes (b,d,g) do not reach the 
on-shell regions when to* is finite. 

For massless quarks, the upper threshold of process (d) and the lower kinematical limits of (b) and (g) are located 
at Ek- However, phasespace just opens at this point and thus contributions to the on-shell width are not generated. 
This picture will not change substantially when collisional broadening is taken into account. If the (on-shell) quark 
width does not become too large - and this will be found numerically, see Fig. [22] - the thresholds are still rather 
sharp. For (b,g) we have also to consider that the constraints from Pauli blocking are stricter than the kinematical 
ones at moderate \k |. 

The processes (a,b,d) of Fig. behave similar to (b,g) of Fig. When is finite, the lower kinematical thresholds 
of (a) and (d) and the upper kinematical threshold of (b) keep the processes off-shell. Pauli blocking prevents (a) and 
(b) even from entering the off-shell regions near the Fermi energy. When m„ vanishes, the thresholds of (a,b,d) will 
be located at Ek- Due to the same arguments as before, there will be no contributions to the on-shell width on the 
quasiparticle level. Small contributions are possible when collisional broadening is considered. 

Process (c) of Fig. [7] is the only process including a bound qq state that may contribute to the (antiquark) on-shell 
width. For this to happen, must be larger than 2to*. The RPA pions of the 0(1/N C ) approach as well as the 
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k Q [GeV] 

FIG. 10: The RPA pion spectral function at a constant three-momentum of 50 MeV. The solid line shows the result of the 
Hartree+RPA approximation in vacuum (m* = 336 MeV), the dashed line the result of the 0(1/N C ) approach at a density of 
0.05po (jri* = 326 MeV). Parameter set I of Tables HI and IIVI has been used. 

physical pions will exceed this limit only in the chirally restored phase. Since the masses of 0(1/N C ) RPA pions and 
Hartree+RPA pions are much closer to each other in the chirally restored phase (we will discuss the numerical results 
for the RPA pion mass in Section IV E[) . such on-shell contributions are no artifact of the present approach. Note that 
the limit k < E^ +kp — /i in Table [TTI is a simplification of k < [i — \(\k\ + k F ) 2 + TO 2 ] 1 / 2 (for \k \ > to* /(m^ — m*)k F ). 

This condition forbids on-shell contributions above a certain value of |fe|. An interesting observation holds for all 
populated on- and off-shell quark states: The contributions to the width - process (b) of Fig. [7] and (c) of Fig. [9]- do 
not include antiquarks. Even when process (c) of Fig. [7] contributes to the antiquark on-shell width, this will have no 
direct influence on the quark states. 

The analysis of the thresholds shows that the on-shell quark width in the chirally broken phase is determined by 
the processes of Fig. GO All those processes do not involve bound qq states (treating the broad RPA sigma peak as 
part of the qq continuum). The thresholds of the processes are independent of the RPA pion mass, cf. Table Hill The 
processes involving bound qq states do not contribute to the collision rates in the on-shell region, even when the pions 
have a reasonable mass. The gap between the on-shell energy and the thresholds of these processes is on the order 
of m w . Considering a realistic pion mass of 140 MeV, this gap is much larger than the typical width of the on-shell 
peaks that we find numerically. 

So far, we have not considered the influence of the pion mass on the magnitude of the various contributions to the 
on-shell width. In the processes of Fig. [9l we do not probe the on-shell region of the RPA pion propagator (i.e., the 
bound qq states). However, the propagator may also depend on the pion mass in its off-shell regions. Thus, it is 
not unlikely that the pion mass has influence on the on-shell quark width even though the width is generated by the 
"correct" processes. 

Let us split up the relevant processes of Fig. [9] into the t-channel processes (a,c,e) and the s-channel process (f). In 
the s-channel process - like in the decays (b,d,g) that do not contribute to the on-shell width - the exchanged meson 
carries a timelike four-momentum that must exceed k 2 > (to* + /i) 2 > (2m*) 2 . In other words, the magnitude of 
the scattering rate is sensitive to the qq continuum of the RPA spectral functions. We show numerical results for the 
qq continuum of the RPA pion spectral function in the Hartree+RPA approximation and in the 0(1/N C ) approach 
in Fig. 1101 The differences in size and shape of the continuum are within reasonable limits - we are comparing a 
leading order to a next-to-leading order scheme here. A significant dependence of the qq continuum cannot be 
observed. This means that the magnitude of the contributions from process (f) to the on-shell quark width is - like 
the thresholds - independent of the large RPA pion mass (in the chirally broken phase) of our 0(1/ N c ) approach. 
We note that the same is true for the contributions of the decays (b,d,g) to the off-shell quark width. 

In the t-channel processes (a,c,e), the mesons must carry a spacelike four-momentum k 2 < [40l |. Typically, one 
would expect that the range of the meson-exchange interaction in these processes scales with the inverse mass of 
the exchanged meson, i.e. R ~ I/to^. The simple estimate T ~ pov, a ~ 7ri? 2 (where a is the total quark-quark 
cross section and v the relative velocity of two colliding quarks) suggests a quadratic to w dependence of the quark 
width, r ~ 1/rn^. Since the RPA pions of our approach are almost three times heavier than physical pions in the 
chirally broken phase, this means that the contributions from the processes (a,c,e) to the on-shell quark width could 
be suppressed. 

The off-shell properties of our dynamically generated mesons differ from those of regular mesons - we come back 
to this point in Appendix [Cj Hence, we cannot directly transfer the above estimate to our approach. A numerical 



19 



check of the term | A 1 ° t (k 7r )\ 2 in (|49)l . that we do not discuss in detail here, shows that the real suppression effect is 
k\ dependent (we integrate over k^ — r — k in the total collision rates) and ranges - on average - well below one 
order of magnitude. The naive estimate T ~ 1/wJ yields only an upper limit for the suppression and overestimates 
the influence of the large pion mass on the on-shell width. 

This brings us to the end of our investigation of the m^ dependence of the quark width. We have found that our 
approach generates the on-shell width from the "correct" processes, i.e. from the same processes as an approach with 
a more reasonable pion mass - this is important for the structure of the width and its density dependence. However, 
the contributions from several processes are suppressed by the too large RPA pion mass. In summary, this means 
that our approach represents a conservative approximation of the short-range effects in quark matter. A calculation 
with a more realistic pion mass should lead to stronger short-range correlations. Let us clarify once again that these 
considerations affect only our results in the chirally broken phase. Above the phase transition, the RPA pion mass of 
our approach will be close to the pion mass of the Hartree+RPA approximation. Thus, there will be no suppression 
of the quark width in the chirally restored phase. 

For completeness, we will consider the correctional diagrams (c,d) of Fig. [6] A self-consistent inclusion of these 
diagrams will fix the pion mass in next-to-leading-order and shift the limits of the processes in Fig. [Jj closer to the 
Hartree+RPA values in the chirally broken phase. More important than this shift, the contributions from the t- 
channel processes of Fig. [9] will rise. Since the chiral properties of the pions remain disturbed in higher orders, the 
RPA pions will still be heavier than physical pions. The correctional diagrams will also modify Eq. (|4"9")l and thus the 
qq continuum entering Fig. [5] These corrections are, however, only next-to-leading order contributions to the quark 
width and thus suppressed in 1/N C . 

We note that the dependence of the thresholds in Table |H] may lead to differences between calculations with 
massive and with massless pions in the chirally broken phase. In contrast to the density dependent (on-shell) processes 
(a,c,e,f ) of Fig. [91 the decays (a,d) of Fig. [7] occur even in vacuum. Thus, their small contributions to the width in 
the on-shell regions - for = may be larger than those of the actual on-shell processes at low densities. 

When moving close to the on-shell regions, the large off-shell contributions of the decays may also have influence on 
the results, even when the on-shell width does not change too much. The left panel of Fig. [T2"l shows a similar effect: 
In the 0(1/N C ) approach, the contributions from process (c) of Fig. [JJ come rather close to the antiquark on-shell 
region in the chirally broken phase. They generate the small peak in the width at low negative kg. This peak is 
reflected in the spectral function by a significant structure next to the antiquark on-shell peak. 

The processes (a) and (d) of Fig. [JJ are less density dependent than (c). Thus, they generate larger off-shell 
contributions at low densities - they can be identified with the shoulders on the large bumps in the left panel of 
Fig. 1121 When those processes move closer to the on-shell regions - which happens when m w becomes small compared 
to the quark width - the impact will surely be more extreme than that of process (c) in Fig. 1121 It is important to 
stress that the pions must be (almost) massless to observe the just discussed effect. The mass of the physical pions, 
m-a ps 140 MeV, is considerably larger than the typical on-shell quark width in our approach, cf. Fig. [^U 

We come to the end of our quasiparticle analysis of the scattering processes. The investigation of the thresholds 
and the density dependence of the individual processes in Sections IIVBI and IIV CI respectively, will be very helpful 
for the understanding of the results of our full 0(1/N C ) calculations. 

V. NUMERICS AND RESULTS 
A. Details of the calculation 

Using the present 0(1/N C ) approach, we have performed numerical calculations for flavor symmetric quark matter 
in thermodynamical equilibrium. We have investigated a wide range of chemical potentials /i, allowing for quark 
matter both, in the chirally broken and restored phase. The calculations have been restricted to zero temperature. 

In the Hartree+RPA approximation, the NJL coupling G and the cutoff A are fitted to the quark condensate and 
the pion decay constant. If a finite current quark mass mo is desired, it is used to adjust m w . Presently, we do 
not calculate the 1/N C corrected physical pions (Goldstone modes) and have no access to their mass and the pion 
decay constant. Hence, we proceeded in a way analogous to the transition from the Hartree to the Hartree-Fock 
approximation in [l3j to adjust the parameter sets of Table U to our approach: We kept A and mo fixed and lowered 
the (Hartree) coupling G by 22% so that a reasonable value for the quark condensate (|32[) was obtained. Note that 
the too large RPA pion mass cannot be fixed by adjusting mo, cf. (|44IA6|) . The readjusted 0(1/N C ) NJL parameters 
are shown in Table HVl 

As discussed in Section ITVl we identify quark states at negative fco with antiquark states k — (fco, k ) at positive fco. 
Note that, strictly speaking, the quark-antiquark border is not located at fco = but at Refco = fco — Re T, Q ct (ko, k ) = 0. 
Re Sq * shifts the energy scale by less than 5 MeV in the chirally broken and less than 50 MeV in the chirally restored 
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TABLE IV: Parameter sets for the self-consistent 0(1/N C ) approach, based on the Hartree parameter sets in Table HI A is a 
three-momentum cutoff. The couplings were reduced by 22% to readjust the quark condensate, m* and {uu) are the effective 
mass ((50} and the quark condensate (|32|) that are found with those parameters in vacuum. 





mo 


GA 2 


A 


* 

m 


/- \l/3 

(uu) 




[MeV] 




[MeV] 


[MeV] 


[MeV] 


Set 




1.68 


653 


308 


-248 


Set I 


5.5 


1.72 


631 


326 


-244 




phase. In the calculations we use the correct border, e.g., when determining the quark density and the momentum 
distribution. In the discussion, however, we just refer to positive and negative energy states to denote quarks and 
antiquarks, respectively. 

For the discussion of the results (not for the numerics) we also introduce an effective quark mass that is real and 
constant. Like in [2(J HH, we define m* in terms of the real part of the on-shell self energy, cf. Eqs. ()19|28|) . for a 
quark at rest, 

to* = m + Re Ef* m (k = 0) = Rem rot - av (fcS s (0), 0) . (50) 

The numerical results show that the momentum dependence of Hem(k) ranges below 5% for the on-shell states. Hence, 
to* yields a good classification for the states in the vicinity of the on-shell peak at all momenta. 



B. Numerical implementation 

The self-consistency problem of Fig.|4]can be solved iteratively. Starting from a Hartree propagator that is modified 
by adding a small width, we calculate the Dyson-Schwinger equations consecutively (beginning with the lower one) 
until the results converge. In the calculations we keep /i fixed while the quark mass and the density may change. 

Fig. [H] shows how the quark width converges in a calculation in the chirally broken phase. The strong variations in 
the first iterations are related to fluctuations of the quark mass. The quasiparticle relation p ~ (p 2 — to* 2 ) 3 / 2 , cf. (|34|) . 
indicates that small changes of m* have significant impact on the density - and thus the width, cf. Section IIV CI - 
when to* and /i are of comparable size. Self-consistency is reached after 10 iterations. 

Calculations in the chirally restored phase, where m* is much smaller than /i, converge faster (4-5 iterations). Close 
to the phase transition, the gap equation may have several solutions. More iterations can then be necessary to find 
one of the stable solutions. 

Solving the width integrals (|37l4ip is the most time consuming part of the numerical calculation - the integrands 
consist of products of spectral functions with pronounced on-shell peaks. We use the package CUBPACK [13] for this 
task. Since it handles vector integrands we can calculate all Lorentz components of the quark width in one turn (per 
iteration). To support the numerical routines, we use a substitution of the general formp 2 — > z — arctan[p 2 — p 2 s (po)] 
for the quark three-momenta. The slope of arctan(x) is largest at x — 0. Hence, the on-shell peaks are smeared out 
in the coordinate z. 
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The quark self energy and the meson polarizations depend separately on the energy and the modulus of the three- 
momentum. Thus, the Dyson-Schwinger equations must be solved on numerical grids. For the quarks and mesons we 
cover ranges of |fc | < 3.5 GeV, \k | < A and |fco| < 3 GeV, |fc \ < 2A, respectively. Using grids with 300 x 100 (quarks) 
and 300 x 200 (mesons) mesh points we achieve a resolution of dko ss 20MeV and d\k | w 6MeV. This choice is 
justified by the results - the structure of E rct and II[ ct is smooth with respect to the mesh size. The spectral functions 
(|20|25|) are not discretized on the grid since the on-shell peaks would be lost. 

Note that the limits for \k | are fixed by the cutoff scheme, cf. Section Ul CI In comparison to these limits, the energy 
limits are chosen rather high. Re E ret and Re n ret are calculated from dispersion integrals that range to infinitely high 
energies. Hence, r SiAI and r CTj7r must vanish at the grid borders. We have to go beyond the thresholds of Section llVl 
here since the quark width generates (small) contributions above the quasiparticle thresholds. The quark spectral 
function at such high energies will not sizably influence the widths (|37l4ip . At higher energies, the on-shell peaks are 
located at three-momenta above the cutoff. Only the much smaller off-shell contributions enter the integrals. 



C. Collisional broadening 

We begin with a qualitative overview. Fig. 1121 shows the quark spectral function Ao, the corresponding width To, 
and Re E et in the chirally broken and restored phases as cuts at a constant three-momentum. The 70-components 
correspond to the norelativistic limits of A and S rct and show all relevant features. Furthermore, Tq yields the largest 
contribution to the width W IT7|) and the quark density |33|) is determined by Aq. 

Note that we use po = 0.17 fm~ 3 to normalize the quark density. The numerical value of po corresponds to 
the nucleon density in normal nuclear matter p nm . Nonetheless, p is here a quark density and not a nucleon density. 
Assuming that each nucleon consists of three constituent quarks, quark matter with a density of po should be compared 
to nuclear matter at p nm /3. 

The spectral functions in Fig. [T2]are dominated by the on-shell peaks. In the upper left panel of Fig. [T21 we find 
the antiquark peak at ko ps — 290 Me V and the quark peak at fco ~ 290 MeV, cf. p^|) . In the right panel, the peaks 
are located at lower energies since the mass of the quarks drops significantly in the chirally restored phase. The 
peak-like structure right next to the antiquark peak (fco ~ —180 MeV) in the left panel (it turns into the shoulder of 
the antiquark peak in the right panel) is generated by a pronounced contribution to the quark width - see the lower 
panel. In contrast to the on-shell peaks, it does not correspond to a pole of the quark propagator. We will come back 
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k Q [GeV] k [GeV] 

FIG. 13: Cuts of the quark spectral function Ao and the width To in the chirally broken phase (p = po) at several three- 
momenta, using parameter sets (dashed lines, [s, = 308.5 MeV, m* = 263 MeV, /cf = 188 MeV) and I (solid lines, /i = 325 MeV, 
m* = 284 MeV, k F = 186 MeV). Note the different scale in k a as compared to Fig. OH 

to this structure in Section [V Dl 

Most of the strength of the spectral functions in Fig. [TJ] is located close to the on-shell peaks, note the log scale 
for Aq. The broad underground of off-shell states is generated by the quark width Tq. A detailed discussion of the 
structure of the width - that also explains the off-shell structure of the spectral functions - can be found below. Here, 
we only note that Tq is strongly energy dependent, far off-shell it rises up to 400 MeV. When the density increases, 
To becomes larger at low |fco| while it drops at large |fco|. As we will see, this is in agreement with the estimates of 
Section [IVS1 

The impact of the cutoff at large |&o| can be clearly observed in Fig.[TJJ To and Ao drop rapidly and will eventually 
vanish. It can be checked that Re£g et approaches S^ ff , cf. Appendix [Cl The influence of the cutoff - that has been 
introduced for technical reasons - is not entirely unphysical. The combination of a constant coupling with a cutoff 
roughly resembles the running coupling of QCD ,1!] : Quarks with large momenta are free and do not participate in 
interactions. Hence, the cutoff can be understood as a crude approximation of asymptotic freedom. In Fig. 1121 the 
width breaks down for the asymptotically free states at large |A»o | ■ 
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k Q [GeV] k [GeV] 

FIG. 14: Cuts of the quark spectral function Ao and the width To in the chirally restored phase (p = 4.5po) at several three- 
momenta, using parameter sets (dashed lines, /j, — 318 MeV, m* = 0, k F = 312 MeV) and I (solid lines, n — 332 MeV, 
m* = 100 MeV, k F = 315 MeV). 



Figs. [13] and [14] show Ao and Tq in more detail for \k \ < E A = k% s (A) g VA 2 + m* 2 . This is the region where 
the on-shell peaks are not suppressed by the cutoff in the integrands of (|37l4ip , cf. Section III CI Cuts for several 
three-momenta are shown, at the same densities as in Fig. [TJ] 

A noticeable feature of the width and the spectral function, are the zeros at fco = Due to Pauli blocking, the 
states at the Fermi energy are stable for T = 0, cf. Section llll Bl Close to the Fermi momentum, the quark peak of Ao 
turns almost into a sharp quasiparticle peak - see the 200 MeV cut of Fig.[T3](/cF = 188 MeV) and the 300 MeV cut 
of Fig. [TJ] (hp — 312 MeV). Note that Aq an d To would not vanish at k = [i at finite temperatures. In vacuum, Ao 
and To would be symmetric in fco while Re Eq * would be antisymmetric. This symmetry is induced by the invariance 
of the vacuum ground state under charge conjugation, see Appendix E in [48| . Figs. 1 1211 141 show that a finite quark 
density breaks the symmetry between quarks and antiquarks - in particular for \ko\ < A. 

The results for parameter sets and I differ only weakly in the chirally broken phase. G and A are of similar size, 
the value of uiq is negligible when m* is large. The differences are slightly larger in the chirally restored phase since 
the effective masses differ significantly: A finite mo breaks chiral symmetry explicitly. Hence, the quarks have a mass 
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TABLE V: The thresholds of Tables [TT] and [TTTJ adjusted to the results for \k \ = 100 MeV and parameter set I in Figs. [121111 The 
left column corresponds to fi = 325 MeV (m* = 284 MeV, m n = 444 MeV), the right column to fi = 332 MeV (m* = 100 MeV, 
m,r = 416 MeV). See the text for details. 

Process Thresholds [MeV] 



ch. broken phase ch. restored phase 







(a) 


769 


< 


fco 


< 


1547 


748 


< 


fco 


< 


1480 


Fig 




(b),(c) 


127 


< 


ko 


< 


189 


136 


< 


k 


< 


331 






(d) 


735 


< 


ko 


< 


1547 


526 


< 


ko 


< 


1480 






(a) 


325 


< 


ko 


< 


1100 


332 


< 


k 


< 


1178 






(b) 


934 


< 


ko 


< 


2076 


764 


< 


ko 


< 


1917 


Fig 
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(c),(d) 


99 


< 


ko 


< 


325 


-73 


< 


ko 


< 


332 






(e),(f) 


284 


< 


ko 


< 


1100 


100 


< 


ko 


< 


1178 






(g) 


893 


< 


ko 


< 


2076 


532 


< 




< 


1917 



of 100 MeV for parameter set I (see also Fig. \T§\\ while they become massless for set 0. Consequently, the on-shell 
peaks are shifted, cf. (|19p . Since we compare results at the same density - not at the same chemical potential - the 
zeros at ko = /i are also not exactly in the same position, cf. (|34[) . 

Fig. [TS] shows a comparison between the present work and our earlier approach in [Tfj| . In the loop-expansion of 
[l(]| we have calculated the quark width from the Born diagrams shown in Fig. [2)Jc,d). Since we found rather small 
widths in jlC| |. it was justified to use the Hartree parameter set of Table [J there. Note that the chemical potentials 
are not equal since the (density dependent) real part of the quark self energy - in particular Re £g 6t - were ignored 
in El. 

Qualitatively, there is some agreement between the results in Fig. [T5J We did not introduce RPA mesons in [10| , 
processes involving bound qq states did not contribute to the width there. Hence, the bump at negative energies is 
only found in the 0(1/N C ) results (see the discussion below for details). Quantitatively, the results differ by one order 
of magnitude. The simple expansion in terms of the coupling constant has clearly missed large contributions to the 
short-range correlations by ignoring the terms that also are of next-to-leading order in 1 /N c but of higher orders in 
the coupling. 



D. Correspondence to scattering and decay processes 

To understand the full structure of Ao in Figs. [T2TfT4l we have to analyze the structure of Tq. Recall our discussion 
from Section HVl - all processes that contribute to the width are shown in Figs. [7] and [5] We will refer to them as 
EH a )j etc. in the following. The energy thresholds of the processes are given - in terms of m* and - in Tables 
III1 and IIIIl In Table [V] we list the numerical values that are found by inserting the masses of the present approach. 
The real parts of E^ Gt have not been taken into account. Hence, the thresholds are slightly shifted with respect to the 
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curves in Figs. fl2lfT4l Note that the strict quasiparticle thresholds are softened by the finite quark width. 

The processes |H{b) and (g) generate the huge bumps of T at large |fc | m Fig. [T3J Such decays can occur at zero 
density since no partner from the medium is required. The bumps decrease at higher p since the outgoing quarks are 
subject to Pauli blocking. The smaller bump at low negative fco is generated by[7Jc). Since the RPA pion is rather 
heavy, the thresholds for[7fb,c) are located at negative fco and[7jb) is suppressed. 

In the chirally broken phase, the decays [7ja,d) can be identified with the shoulders on the inner slopes of the 
big bumps in Fig. 1121 [T^a) remains visible in the chirally restored phase while [7£d) is hidden below the increased 
contributions from[5Je,f). [DJe,f) contribute to Tq between the large and the small bump in the antiquark sector. The 
contributions from [9ja) and 02c, d) are located right above and below the Fermi energy, respectively. Since [9ja,c-f ) 
depend on quarks from the medium, their influence on r increases substantially in the chirally restored phase. 

Note that [Tfc) generates a structure in Aq that looks like a second antiquark peak in the chirally broken phase at 
fco ~ — 180MeV, cf. Fig. [T2J This is an artifact of the present approach. For a more reasonable in the chirally 
broken phase, the contributions from [Tie) would be located further off-shell. Ao at positive energies is free from such 
artifacts. A process comparable to Etc) ,i.e. kq — > n, would require an antiquark from the medium. Fig. [14] shows 
that process W[c) has considerable influence on the antiquark on-shell width in the chirally restored phase. There, 
however, the Hartree+RPA pion acquires a considerable mass, too. Thus, the contributions of [TJc) should be located 
approximately in the right position. 

The quarks of the medium populate all states at positive energies below the chemical potential. In the cuts for 
\k\ > kp in Figs. [T51 and [T4l the quark on-shell peak is located above the chemical potential. Nonetheless, the 
spectral functions do not entirely vanish in the range < fco < /i. Populated off-shell states exist due to the width 
generated by[9^c). Those states are important for the momentum distribution of the medium that will be investigated 
in Section IVHl 

We close the discussion on the structure of To and Ao with an important observation: The numerical results confirm 
our quasiparticle considerations of Section [TV] The width in the range |fco| < E\ (« 630 — 700 MeV) - where the 
on-shell peaks are not suppressed - is mostly generated by the processes of Fig. [9] Those processes do not involve 
bound qq states. It is interesting to note that the contributions of the collisional processes to the overall width in the 
chirally broken phase are small in comparison to the huge off-shell contributions from decays. This can be readily 
explained by the strong density dependence of the collisions with partners from the medium. In addition, we have 
found in Section HV Dl that the contributions from the t-channel (a,c,e) processes of Fig. |H] depend - in contrast to the 
decays - on the RPA pion mass. Thus, we may underestimate these processes below the chiral phase transition. 



E. RPA mesons in the 0(l/N c ) approach 

The RPA mesons that we discuss here cannot be identified with the physical a and tt - the \/N c corrections of [23|, 
cf. Section IIII El are presently not calculated. Fig. [T|)] shows the spectral functions and widths of the RPA mesons 
in the chirally broken phase. Since parameter set is used (chiral limit), the physical pion should be massless. In 
the 0(1/N C ) calculation, however, we find m a — 603 MeV and m w = 424 MeV. As expected, the RPA pions are no 
Goldstone bosons. 

The widths in Fig. [Tulare dominated by the qq continuum. The quasiparticle threshold for qq decays is fco = fi+m*, 
above 2Ea = 2 (A 2 + m* 2 ) 1 / 2 the decays are suppressed by the cutoff. We find small but finite widths beyond these 
limits in the 0(1/N C ) calculation. The magnitude of the qq continuum drops for larger \k |. For |fc | > 2A, at least 
one quark propagator in (|39p must have a momentum above the cutoff. Consequently, the integral becomes zero and 
the qq continuum vanishes. 

At spacelike four-momenta, Landau damping (qir — > q, qa — > q) [43[ generates another structure. Since the widths 
are given by T lTj7r (k) = — Imn™^(fc)/fco (j30|) . they can become rather large at small fco- However, they do not diverge 
for fc — > 0. 

The continuum parts of the spectral functions £ CT)3r in Fig. [Tulare in good agreement with the established mean field 
approaches. They should generate reasonable contributions to the quark width. The on-shell peak of the RPA pion 
- generated by bound qq states - is located below the qq continuum. Thus it is rather sharp in comparison to the 
sigma peak. 

In Fig.[T7]we take a closer look at the real parts of the denominators of A™^.. The zeros of 1 + 2GRe T1™1 determine 
TO CT 7r . At small fco, the 0(1/N C ) pion curves are shifted by ~ 0.2 in comparison to the Hartree+RPA result. This is 
in agreement with the estimate from Section IIII El About 50% of the shift can be attributed to the Fock self energy, 
cf. (|A6|) . Dispersive contributions to Re n ret generate the remaining 50%. increases from 142 MeV (Hartree+RPA) 
to 460 MeV. Note that the RPA pion is not only heavier than the physical pion but also broader. The on-shell peak 
is shifted closer to the qq continuum, where the width is larger. 
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FIG. 16: Cuts of the RPA meson spectral functions and widths at constant three-momenta. The calculation was performed in 
the chirally broken phase at \i — 308.5 MeV (p = po, m* — 263 MeV), using parameter set 0. 
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FIG. 17: The real parts of the denominators of at \k \ = 50 MeV in the chirally broken phase. The solid lines show 

Hartree+RPA results at zero density, the dashed and the dotted lines results from 0(1/N C ) calculations at p = 0.05po and 
p = po, respectively. Parameter set I from Tables HI and IIVI has been used. 



In the case of the RPA sigma, the shift of the real parts is smaller. Furthermore, the zero is located in a region 
where the shift has - due to the slope of the curve - less influence on m a . The 0(1/N C ) effects are only on the order 
of 10%. Note that Fig. [T7] also shows that m a has a stronger density dependence than m,. 



F. Chiral phase transition 



To examine the influence of the short-range correlations on the chiral phase transition, we investigate the effective 
quark mass m* ([50]) as a function of the chemical potential. Before discussing the 0(1/N C ) results, we first take a 
closer look at the Hartree(-Fock) results in Figs. [T] (solid line) andfTSl 

In the shaded region of Fig. [H the gap equation (cf. Section Hi C| has three solutions m*(fj.) for a given value ofu. 
The (inverted) 'S' shape of the curve is characteristic for the order parameter of a first-order phase transition [T{|. 
The upper and the lower solution (chirally broken and restored phase) correspond to minima of the thermodynamical 
potential [l3l |50] | and thus are stable. The intermediate solution corresponds to a maximum of the potential and 
is metastable. Note that m* remains finite in the chirally restored phase since a small mo breaks chiral symmetry 
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FIG. 18: The effective quark and RPA meson masses (left) and the quark density (right) in the Hartree+RPA approximation, 
using parameter set I from Table [U A first-order chiral phase transition occurs at p = 355 MeV. (Data taken from [l3L 14^.) 
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FIG. 19: The effective quark and RPA meson masses (left) and the quark density (right) in the 0(1/N C ) approach, using 
parameter set I. In the range 326 MeV < fi < 328 MeV a first-order phase transition occurs. 



explicitly. 

The energetically favored solution for to* corresponds to the lowest minimum of the thermodynamical potential. 
Depending on fi, this can be either the upper or the lower solution. The phase transition is located where both minima 
have the same value. At this point, to* drops discontinuously while p, cf. (|34[) . increases. The intermediate range of 
(metastable) masses and densities is not realized. The quark density p — 3po - corresponding to nuclear matter at 
normal density - lies within this gap in Fig. 1181 This is a known shortcoming of the NJL model on the quasiparticle 
level 0. 

Fig- US shows the effective quark mass that is found in the 0(1/N C ) calculation, using parameter set I. In contrast 
to the simple estimate in Fig. [T] (dashed line) that ignores the energy and density dependence of the short-range 
correlations, the phase transition is still of first order. The 'S' shape in Fig. [TH] is interrupted since the iterative 
calculations will not converge for the metastable solutions in the vicinity of the phase transition. 

The short-range correlations are strongly density dependent since they are generated by interactions with the 
medium. Only when p - and thus p - increases, deviations from the mean field models can be observed in Fig. [T9j 
The right panel of Fig. [19] shows the strong p dependence of the quark density. Recall the quasiparticle relation 
p ~ (p 2 — m* 2 ) 3 / 2 , cf. (134[) . The explicit cubic p dependence of p is enhanced where m*(p) changes rapidly, i.e., near 
the phase transition. 

The region of chemical potentials for which several solutions for to* exist shrinks from 4 MeV (Hartree) to only 2 MeV 
in the 0(1/N C ) approach (we will not determine the exact location of the phase transition from the thermodynamical 
potential). For the gaps between the upper and the lower solutions for to* and p we find Am* = 111 — 116 MeV and 
Ap = 2.1 — 2.2 p - The Hartree approximation, on the other hand, yields Am* = 188 MeV, Ap — 4.4po [HI- 

We conclude that the short-range correlations are not strong enough to turn the first-order phase transition into a 
smooth crossover at zero temperature. However, the gaps at the phase transition have decreased by 40% — 50%. The 
quark density of 3po is still located within the gap. Note that the smaller to* in the chirally broken phase and the 
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FIG. 20: The quark condensate as a function of the density. The gap in the curves corresponds to the range of densities that 
is skipped at the phase transition, cf. Fig. 1191 

0(1/N C ) term E^ ff , cf. (|C5IC6|) . shift the phase transition to a lower \i. 

Figs . [TBI and [T9l also show the RPA meson masses. The RPA pion of the G(l/N c ) approach has a considerable mass 
in the chirally broken phase. The RPA sigma mass is closer to the Hartree+RPA result. At the phase transition, mo- 
does not drop as far as in the mean field calculation, m n drops instead of increasing, thus moving towards the mean 
field result. At higher fj,, m a and m n resemble the Hartree+RPA result qualitatively and quantitatively. m a and m w 
converge at large fi (chiral partners) - restored chiral symmetry is still a good symmetry in the 0(1/N C ) approach. 

We do not show the results of a calculation using parameter set explicitly here. In the chiral limit (mo = 0), 
the results for the chirally broken phase - including - are very similar to those in Fig. 1191 G and A are almost 
equal in both parameter sets, mo is irrelevant when m* is large. Since E^ 6 ', and thus m* , vanishes above the phase 
transition in the chiral limit - chiral symmetry is fully restored - the 'S' shape of the m* curve in Fig. 1191 turns into 
a '2' shape. Note that the solution m* = of the gap equation exists for all /i in the chiral limit. At large fi, where 
m* becomes very small in Fig. 1191 the results for the two parameter sets are close to each other again. 

Fig. [2TJ1 shows the quark condensate (uu). Like m*, (uu) is large at low densities. The results in the chirally broken 
phase are almost equal for both parameter sets. When parameter set is used, the quark condensate drops to zero 
above the phase transition. It drops but remains finite when chiral symmetry is explicitly broken. 

G. On-shell self energy 

For a quantitative investigation of the collisional self energy we will restrict ourselves to the on-shell states (|19[) . 
Fig. shows that most of the strength of the spectral function remains in the vicinity of the on-shell peaks. The 
off-shell structure of the self energy has only limited influence on the properties of the medium. 

Our calculations show that the energy and momentum dependence of Rem os - the real part of m ret (&;) (fl5|) on the 
mass shell - is weak. Due to the large constant Hartree self energy, the effective masses at k = and \k \ = A differ by 
only 4 — 5%. The Fock self energy E^ = E H / (2NfN c ) that is part of the meson exchange provides another constant 
contribution. Adding up mo + E H + Ef , we find that 80 — 90% of Rem os are determined on the mean field level. Only 
the remaining 10 — 20% are generated by short-range correlations. Recall that the coupling - and thus E H (|3Tj) - has 
been lowered by 22% in the 0(1/N C ) approach. Hence, the collisional self energy just replaces the missing mean field 
contributions without generating new effects. 

The 7^-components of ReE^* shift the peaks of the spectral function, cf. (fT6| . We will not discuss ReEJf* s here 

since it is small in comparison to the other components and to \k\. Re Eo C * s is shown in the upper panels of Fig. 1211 

In the chirally broken phase, Re Eo C * s is moderately density and energy dependent. Up to a small shift (Ef ff is small 
at low fx, cf. Appendix [C| it is approximately antisymmetric in ko. It cannot be absorbed into \i like the constant 
Re E ct = Eq of the Hartree-Fock approach. 

Re E C o S becomes larger in the chirally restored phase. The energy dependence increases and the antisymmetry is 
lost. The non-dispersive Eg ff is on the order of 10 — 25 MeV. At low energies, Re E™os has the same sign for quarks 

and antiquarks. Thus, both peaks of the spectral function are shifted to higher energies, see the cuts for small \k \ in 
Fig. [14] Note that the density dependence has grown for the quarks but is still moderate for the antiquarks. 
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FIG. 21: On-shell real parts of Eq * and the on-shell width for the quark and antiquark states at five different densities in 
the chirally broken (0.5po,Po) and restored phase (4po, 5po, 6po), using parameter set I. The curves start at fco s (|fc| = 0) and 
end at fcg s (|fc | = A). Note the different scales. 



The bump at ko as 120 MeV (right panel of Fig. l2"Tj) is generated by process U^c). Above the phase transition, this 
process moves into the antiquark on-shell region. As discussed before, the processes involving bound qq states have 
no influence on the on-shell self energy in the quark sector. 

Re Eq C * s clearly modifies the properties of the medium in comparison to the mean field approaches where Re Sq C o S 
is zero or constant. Shifting the quark peak of A with respect to the (fixed) chemical potential has, e.g., influence on 
to* and p, cf. (|31l33p . Since RcEq C q S increases at large p, the impact is stronger in the chirally restored phase. 

We turn now to the on-shell width. The lower panels of Fig. [5T] show Tq s - the largest Lorentz component of T os - 
for the quark and antiquark states. Like in the cuts of Figs. 11211141 the width vanishes at ko = fx. In the vicinity of 
fi we find a quadratic energy dependence. Below the phase transition, the populated quark states are located close 
to the Fermi energy. Their width - generated by process - remains small. Above the Fermi energy and in the 
antiquark sector, Tg s shows a linear density dependence - the relevant processes EJa) and[9](e) both involve one quark 
from the medium. 

In the chirally restored phase, i.e. for the densities > 4po in Fig. the region of populated states becomes 
larger and Tg s increases up to 10 — 20 MeV. It drops at low ko since the threshold for process |5fc) is approached, 
cf. Section HVBl The widths of the free quark states and of the antiquark states are again larger. The density 
dependence of Tq s weakens drastically at higher p. We can interpret that as a - Pauli blocking and cutoff induced 
- saturation of the short-range effects, cf. Section IIV CI We have observed the same effect for nucleons in nuclear 
matter at nucleon densities of a few times po [46] . 

For a better understanding of the density dependence of the width, we can investigate the average width of the 
populated quark states (r pop ) and of all quark states (rg 11 ) (excluding antiquarks). We define (Fg° p ) as 

<r po P> = Jo A dkk2 fZkp=o dkoTo(k)Ao(k)n F (ko) 

I A dkk2 JZk o =0 dk Q Ao(k)n F (k ) 

using again the three-momentum cutoff A. Note that the denominator is just the density (|33|) . (rg 11 ) is found by 
removing the distribution functions in (|5ip . Since we weight Tq with the spectral function, the on-shell width will 
dominate the average widths. 

Fig. |2"21 shows (Fq° p ) and (r^ 11 ) from calculations with both parameter sets. (Fq° p ) remains smaller than 2 MeV in 
the chirally restored phase. Above the phase transition, we can observe the expected saturation effect - (r pop ) does 
not exceed 10 MeV. (Tg 11 ) grows up to 20 MeV below the phase transition. In the chirally restored phase, however, 
it does not only saturate but starts to decrease. Note that the structure of the curves in Fig. [2U can be readily 
understood from the considerations in Section ITVl We will not discuss it here in more detail. 
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FIG. 22: The average widths (Fq° p ) and (Fq 11 ). The gap in the curves corresponds to the range of densities that is skipped at 
the chiral phase transition, cf. Fig. 1191 

The results for the average widths appear rather small in comparison to Figs. 1131 1 141 and 1211 They are, however, 
comparable to the nuclear matter results in [46|]. At T — and p nm — po, the average width of the populated nucleon 
states is 1.5 MeV. At 3po, a width of 6MeV is found. Of course, nucleon and quark widths should not be compared 
directly. Nonetheless, the quantitative agreement indicates that the short-range correlations are of similar importance 
in both systems. 

In Section IV Fl we have found that the chiral phase transition does not turn into a smooth crossover like in the 
naive calculation of Fig. Q] In the extended Hartree approach, we have used a constant width of 24 MeV. Below the 
phase transition, the realistic value for (Fg° p ) is one order of magnitude smaller. Hence, it is not surprising that we 
still find a first-order phase transition. Recall that we consider a conservative estimate of the short-range effects in 
the chirally broken phase here, cf. Section TlV Dl A calculation with a more realistic pion mass should yield an even 
smaller mass gap or turn the phase transition into a crossover. We note that the nucleon width is strongly temperature 
dependent [461 ]. A similar behavior can be expected for the quark width. This should lead to a stronger influence of 
the short-range correlations on the phase transition at finite T. Such a study is, however, beyond the scope of the 
present work. 

H. Momentum distribution 

To estimate the importance of the off-shell states further away from the on-shell peaks, cf. Figs. fT2lfT4l we investigate 
the momentum distribution of the quarks in the medium. The momentum distribution is closely related to the density 
(|33| . It is given by 

n(\k\) = - I dk A (k)Mko) ■ (52) 

n JRcfc o =0 

The factor 1 /tt has been chosen to normalize the momentum distribution, cf . (|10[) . 

In a quasiparticle approach (|45p . we would find the momentum distribution of a free Fermi gas that turns into a 
step function 0(/cp — |fc|) at T — 0. The finite quark width modifies this simple picture. The off-shell states shift 
some strength away from the peaks in Ao- For \k | < kp this leads to a depletion of the momentum distribution since 
states above the chemical potential are not populated. For |fc| > kp, the off-shell states shift some strength back 
below the chemical potential. Those populated states yield - even at T = - a so-called high-momentum tail in the 
momentum distribution. The depletion and the high-momentum tail can be interpreted as a universal measure for 
the short-range correlations. 

The left panel of Fig. [2"31 shows the momentum distribution at several densities. As expected for an infinite system 
at T = 0, we find a sharp step at Hf Q- Recall that the width in the range < k < fi is generated by process HJc) . 
The density dependence of this process saturates above the phase transition, cf. Fig. [55] (left panel). This behavior 
is reflected in n(\k |). While the high-momentum tails at the lower densities differ significantly in size, they are very 
close to each other at the higher densities. 

At low | A; |, the momentum distribution has values of 0.94 below and 0.89 above the chiral phase transition. The 
discontinuity at Uf has a size of 0.82 — 0.85 for the lower and of 0.72 — 0.75 for the higher densities. These results 



31 




k [GeV] k [GeV] 

FIG. 23: The momentum distribution of quarks in quark matter. The left panel shows 0(1/N C ) results, using parameter set I 
from Table ITVl The solid and the dashed line correspond to quark matter in the chirally broken phase, the other lines to quark 
matter in the chirally restored phase. The right panel shows results from [lOj], using parameter set from Table U and restoring 
chiral symmetry by hand. 

can be readily compared to the well established nuclear matter values, see Q for a compilation. At T = and 
normal nuclear matter density, the depletion of the nucleon momentum distribution is on the order of 15 — 18%. The 
discontinuity at fcp- has a size of 0.7 — 0.75. 

We can conclude that the short-range effects in quark matter in the chirally restored phase are of similar size as in 
normal nuclear matter. In the chirally broken phase, where the quarks reach the same particle number density po as 
the nucleons at normal density, the correlations are slightly smaller. 

It is also interesting to compare the present results to our previous calculations. The right panel of Fig. [23] shows 
the momentum distributions from [To| at two different densities (m* = 0). The high-momentum tail at p = 3po is 
more than one order of magnitude smaller than the 4po result of the 0(1/N C ) approach. A density of 9po is needed 
to generate short-range effects on the same level as in the 0(1/N C ) approach at 0.5po — lpo- This confirms again, 
cf. Section [VDl that the loop-expansion has missed significant contributions to the short-range correlations. 

VI. SUMMARY AND OUTLOOK 

In this work we have explored the content of short-range correlations in cold and dense quark matter within the 
framework of the SU(2) NJL model. Our approach is based on techniques that haven proven to be very successful in 
nuclear matter. Employing a next-to-leading order expansion in 1/N C for the quark self energy, we have constructed 
a fully self-consistent model. By summing up all 0(1/N C ) self energy diagrams, we have dynamically generated RPA 
mesons. Their properties are - in contrast to the Hartree+RPA approach - self-consistently fed back into the quark 
self energy. To preserve analyticity, the reals parts of the quark self energy and the RPA polarizations have been 
calculated from dispersion integrals - non-dispersive terms were identified explicitly. 

The RPA pions are not the Goldstone modes of our approach. Presently, a next-to-leading order approach to the 
NJL model that is fully self-consistent and yields Goldstone pions does not exist. An analysis of the structure of the 
quark width has shown that the too large pion mass does not influence which processes contribute to the on-shell 
width. However, the spacelike off-shell components of the RPA pion propagator, that are important for the generation 
of the on-shell quark width, depend on the RPA pion mass. Thus, our results for the chirally broken phase represent 
a conservative estimate of the short-range effects. 

In calculations at finite chemical potentials and zero temperature, we have investigated the (off-shell) structure of 
the quark width, the size and density dependence of the short-range effects, and their influence on the properties of 
the medium. The short-range correlations are not strong enough to turn the first-order phase transition into a smooth 
crossover at zero temperature. Nonetheless, the mass and density gaps at the chiral phase transition have dropped 
by 40 — 50%. Above the phase transition, the correlations quickly saturate due to Pauli blocking and the NJL cutoff. 

We have found quark widths that are similar in shape but one order of magnitude larger than in the loop-expansion 
of [l(j. The larger results of the present work arise from a consistent expansion of the quark self energy in 1/N C . 
The comparison of the quark matter results to nuclear matter calculations shows that the short-range correlations 
are of similar magnitude in both systems. A complete description of nuclear systems is only possible when the short- 
range effects are taken into account. Therefore, the influence of short-range correlations on the quarks (and the RPA 
mesons) should also not be ignored. 
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Let us consider the QCD phase diagram to stress this point: Phenomenological models are used to determine the 
properties of the chiral phase transition. Many models suggest that a first-order phase transition turns into a crossover 
at a critical point. As we have seen in this work, short-range correlations change the properties of the phase transition 
with respect to quasiparticle calculations. Nuclear matter calculations indicate that the short-range correlations will 
be even larger at finite temperatures [461 ] . Hence, the mean field models will most likely predict a wrong position for 
the critical point. 

The present 0(l/N c ) approach can be improved in many ways. The most interesting but also most complicated 
extension would be to fix the chiral properties of the RPA pions. A first step into that direction would be to include 
the next-to-leading order corrections to the RPA polarizations - or to generate higher order vertex corrections via 
Bethe-Salpeter equations - in the self-consistent calculation. Such an approach would allow us to take a closer look 
at pions with more reasonable masses and at the chiral theorems in the presence of short-range correlations. Our 
investigations from Chapter IIVI indicate that the on-shell width of the quarks could also be larger than determined 
here. 

Most likely, the short-range effects in quark matter will - like those in nuclear matter - increase significantly 
at higher temperature. This should be checked in an explicit calculation. We have already discussed the possible 
consequences for the structure of the chiral phase transition and the critical point. 

We have ignored the phenomenon of color superconductivity flU [ItJ in the present approach. Diquark states can 
be dynamically generated in the same fashion as the RPA mesons. A diquark condensate may have influence on the 
short-range effects in the chirally restored phase [5l| . The short-range effects will also have some influence on the 
diquark condensate. For the investigation of systems like neutron stars and hypernuclear matter, it would also be 
interesting to extend our model to asymmetric quark matter and flavor SU(3). Similar approaches exist already for 
nuclear matter. 
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APPENDIX A: RPA ON THE HARTREE(-FOCK) LEVEL 



For the present work it is instructive to review certain aspects of the quasiparticle RPA. An investigation at 
T = fi = (chirally broken phase) provides insights about the Goldstone boson character of the RPA pions and the 
proper use of dispersion integrals. 

It has been shown in that the time-ordered RPA pion polarization (|B2[) of the Hartree(-Fock)+RPA approach 
at T = /i = has the form 

-iHl(k 2 )=-W f N c / f P - 1 >2 +2JV / JV c fc 2 I(fc 2 ) (Al) 
J (27r) 4 p z — m 

where to* = to + S H (+ ) and 

I(fc2) = / (2tt) 4 [(p + ifc) 2 - to* 2 ] [(p - \kf - m* 2 ] ' (A2) 

Note that the imaginary part of I corresponds to the real part of IT . The time-ordered quark propagator is given by 
S" 1 = # - m*. Using Eq. (01]), 



we can identify the denominator of the RPA pion propagator A£ ([4]) with 

l + 2Gn^(fc 2 ) = m *~ S +^GA^ e fc 2 I(fc 2 ). (A4) 
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If the pions are Goldstone modes, their mass should vanish in the chiral limit (and be small for finite mo) when 
chiral symmetry is spontaneously broken. should then have a pole at k 2 — 0. The second term on the rhs. of (|A4[) 
is surely zero for k 2 = 0. For the first term we find in the Hartree approximation, to* = mo + S H , 

m* m Q + E H ' y ' 

This expression is zero in the chiral limit. For finite too we have too/ E H ~ 0.016. Hence, the Hartree+RPA pions 
behave like Goldstone bosons. 

The Hartree-Fock approximation is not a clean expansion in 1/N C . However, we can use it to estimate the next- 
to-leading order effects. Using to* = to + E H + and S H = 2N f N c the first term of flU} becomes 



to* - E H 1 



E + Eg 2N f N c + 1 



0.08. (A6) 



This term is now significantly larger than too/ E h . Thus, the Hartree-Fock+RPA pions will not behave like Goldstone 
bosons - even in the chiral limit they will have a finite mass. Fig. [5] shows 1 + 2GRen[ ct from a Hartree+RPA and a 
Hartree-Fock+RPA calculation. The zeros determine the pole of the propagator. We can see that even small shifts 
have significant influence on the pion mass. The shift of 0.08 raises by 150 MeV. The 0(1/N C ) approach includes 
the Fock self energy. If no cancellation with other 0(1/ N c ) terms occurs, the 0(1/ N c ) RPA pions will behave even 
less like Goldstone bosons. 

The denominator of is found in the same way as (|A4|) , cf . [l3| , 

1 + 2G n= (fc 2 ) = m * m ^ + HGNfNc (k 2 - 4to* 2 ) I(fc 2 ) . (A7) 

The first term is of minor importance here. As long as it becomes not too large, the poles of A^. et are determined 
approximately by the zero of the second term, i.e., m a w 2m*. 

We turn now to the second term on the rhs. of (|Al|) . The real and the imaginary part of l(k 2 ) can be calculated 
explicitly using either the residue theorem or Cutkosky rules [l[ when a three-momentum cutoff is used to regularize 
the integral. Such a calculation shows that Rel(fc 2 ) sgn(fco) - the main contribution to Imn^, cf. (|B1IA1[) - and 
Iml(fc 2 ) are related by a dispersion integral without additional constant terms, 



~T / ^0 2 _ ,2 

7T Jo % K 



Imr ct (fc , k ) = V dq 2 2 v 7' 1 , ( A8) 



where Rel ret (/fc , k ) = Rel(k 2 ) sgn(k ) and Iml ret (fc) = Iml(fc 2 ), cf. flBTJ. 

When KeU r ( ^ t 7r (k 2 ) shall be calculated from a dispersion integral, it must be taken into account that the ko de- 
pendence of Imll^ is not given by Rel rct alone but by (k 2 — Am* 2 )Rel ret (ko,k 2 ) and fc 2 ReI ret (fc , k 2 ), respectively. 
Calculating the dispersion integral over Rel ret and multiplying the result with a factor k 2 - i.e., inserting the disper- 
sive result for Iml into Re n™^. - is not equivalent to calculating a dispersion integral over fc 2 ReI rct - Imll^. The 
difference is a ko independent term, 

fc 2 lml(fc ,fc) = --k 2 V r ^o 2 ^ 1 ^, ^ (A9) 



7T 



1o ^0 



n Jo t io % ~ K 

i r , 2R Tret , i y r, 2 ininf ( g0 , fc) 

/ dq Rel (g ,fc)- V / dq g -j- 



7T 



const, in ko 



and likewise for the a case. 

Using the dispersion integral over Rel ret and multiplying with k 2 afterwards, is equivalent to the direct calculation 
of Re H™^ . Hence, the ko independent term in the last line of (IA9|) must be calculated explicitly when Re is 
calculated from a dispersion integral over Im U 1 ^ , 



!<n r^ ImI C(<?o,AO 2N f N c f \ 2 ^ TIO w ^ S H 



RelC(*) = ^y o *t P,/ ^I-\ qo ,k)-^. (AlO) 



34 



The first term on the rhs. is the regular dispersion integral. The third term is the constant shift from (|A3[) . A 
decomposition similar to (|A10[> is found for the 0(1/N C ) approach, cf. Appendix [El Note that Il n in (|B22p corresponds 
to the third term of (|A10|) and Ilf in (|B22|) corresponds to the sum of the first and the second term of (jAlOjl . 



APPENDIX B: REAL PARTS OF THE RPA MESON POLARIZATIONS 

In the following we will examine the complex polarizations to identify the non-dispersive part of Re Ilf. Likewise, 
we will investigate the quark self energy in Appendix [C] The Feynman rules of the real-time formalism provide no 
direct access to the retarded self energy and polarizations [3l|, [33| . Hence, we will start from the time-ordered E c and 
H7. They are related to E rot and Ilf in a simple way [33]: 



ReS c (fc) = Re£ rct (fc), 
Rell^fc) = ReIlf (fc), 



Im£ c (£:) = [1 - 2n F (fc )]Im S rot (/c) , 
Imnf(fc) = [1 + 2n B (fco)]Imn[ ct O) 



(Bl) 



with the zero temperature limits 1 — 2n F (fco) — ► sgn(fco — n) and 1 + 2ne(fco) ~~ * sgn(fco). The imaginary parts of S r 
and Ilf that are found in this approach are of course equivalent to the widths that we have found earlier (|37l44p . 



1. Decomposition of the RPA polarizations 



In this section, we will frequently use the relations between the different kinds of propagators (S c ,S ret ,S < ) from 
Section III Dl without further reference. For the rather technical manipulations we introduce the shorthand notation 
p± = p± \k and, cf. (|T5)> . 



n F = n F (p ) , 111 
The time-ordered polarizations of the RPA mesons are given by 



~_rct,av 
± 



m + K Bt ' aV (P±) 



-»D?(*) = - 



d 4 p 



Tr 



r, s c + r, s< 



(B2) 



We can replace the time-ordered propagators by S c (fc) = S rct (fc) + S < (fc), 



n?(fc) = -i 



d 4 P 

(27T) 4 



Tr 



r, sf r, s r _ et +r, sf r, sf +r, sf r, s rct +r t si r, sf 



(B3) 



The first term of the integrand is zero, cf. (HU). Using S rct (fc) = ReS rot (fc) - %A{k) and S<(fc) = iA{k)n F (k ), the 
polarizations can be split up into their real and imaginary parts. With the help of (|B1I) and the relation np(j>+) + 
n F (p_) - 2n F (p+)n F (p_) = [n F (p_) - n F (p+)][l + 2n B (fc)] we find 



nf (fc) 



d 4 p 



Tr 



f 1 s r + et r ; sf +f 1 sf r, s-f v 



(B4) 



Replacing the non-ordered propagators in (|B3|B4j) by S± = [S^ — Sf ]n F yields expressions for n c ' ret that contain 
only retarded and advanced propagators, 



nf (fc) = 



d 4 P 

(2tt)* 



d 4 p 



Tr 



Tr 



-r, sf r, sf et n ; + r, s* v r, s* v n+ 



-f, Sf T, S!l v {n F -n+} 



f , Sf S r _ et (l - n+)(l - n F ) + f , S; v T ; S!l v n+n F 



(B5) 



(B6) 



+f , Sf F, S!_ v (l - n+)n F + f , S* v T Z S r _ et n+(l - n F ) 



35 



These expressions are a good starting point for our search for the non-dispersive part of ReII| ot (fc). The structure of 
gret.av resem ^i es _ ev en in an interacting medium - that of free propagators, cf. (|14[) . Thus, an investigation similar 
to the Hartree+RPA approach at T = /i = 0, see Appendix lAl and [1 31 ] for details, is possible. 

Instead of inserting the explicit forms of S ret ' av (TH)) into (|B5|) or (|B6[) right away, we will investigate the general 
expression TrT; S+ Ti S_ first. In this trace, the indices ± do not only mark the four-momentum but implicitly also 
the kind (ret, av) of the propagators and the self energies contained in them. Later we will derive expressions for 
specific combinations of retarded and advanced propagators from this general form. Inserting (| 14[) into the general 
trace, we get 

TrIV S+ IV S_ = 4N f N c _ 2 + , 2 . (B7) 

(pi — m + )(p_ — m_) 

The upper sign refers to the a and the lower sign to the tt case, respectively. T n is used as a shorthand notation for 
T ± (isospin symmetric matter). Note that the numerator would take the simpler form p 2 — jk 2 ± m* 2 in the Hartree 

approximation where Sg Ct,av (fc) = E mf and Ep e ^ av = 0. 



We can rewrite the denominator of (|B7|) in terms of partial fractions 
Tr IV )7r S+ r CTi7r S_ 
- 2N f N c 



1 



pi — mi pi — m_ 



J (^-t?i^(p 2 _-m 2 _) V ; 



The first term on the rhs. corresponds to a constant contribution to ReII[ ct . It can generate a pion mass and may 
spoil the Goldstone boson character of the RPA pions. The second term is complex. Its real part can in principle be 
calculated from a dispersion integral. k independent shifts that might arise from this term will be discussed later. 

Inserting the first term on the rhs. of (|B8j) into the retarded polarizations (|B5|) yields for the non-dispersive, constant 
part of the meson polarizations 



n n = -UN f N c I -^(-2i)Im 1 n F (p ) . (B!) ) 

\ Z7T ) Piet ~ m ret 



To obtain this expression, we have used p* ot = p av , fh* et — m av and have performed the substitutions p — > p ± \k. 
The rhs. of (|B9|) resembles the Hartree self energy (|3"Tj) . Introducing the constant term E^onst we nn( i 



E^-ieGiV/iV^mo + Seonst.) / WTi 1 ^ 1 ~2 MPo) 

J Prct m rct 



16GN f N c I -JLlm S J > ™ nSt - MPo) ■ (BIO) 



(2tt)4 p 2 



ct 



The value of S* onst can in principle be freely choosen to shift strength between the two terms. One choice is to collect 
all constant contributions to y,° s ' ret in Sconst.- Inserting (|B10j) into (|B9|I yields 



, S mf , 8N f N c f d*p Sf (p)-Konst. , , min 
H n = ~^7T- I ^ + ~ I ^ / 7^34 Im — 3 ^2 Mpo) ■ (Bll) 



const, , 

2Gm +E" onst . ' m + Konst. J ( 27r ) 4 "' Pre 



ret 



The first term on the rhs. corresponds to the result of the Hartree(-Fock) approximation (IA4j) when S* onst is set 
to S H (+S F ). The second term results from using full propagators. In the Hartree(-Fock) approximation, where 
E^ et (fc) = S H (+S F ), this term will vanish. 



3G 



We turn now to the second part of 
the time-ordered polarizations 



n^(fc) = 2iN f N c 



2. The dispersive part of the RPA polarizations 

]B8[) . Inserting this expression into Eq. (|B6|) we find for the dispersive part of 



d 4 p 

Wf 



■p r _ et ) 2 - (m^ et ±m r _ et ) 



(P4 



m+ t2 )(p r _ et2 



g-(l-n+)(l-iiF) 



(m+ ± m a 







- m av2 ) 


(pf - 


pay) 2 „ ^ret 


± 77li V ) 2 


(pf 2 


- mf 2 )^ 2 


- TO^ 2 ) 


(p+- 


p r _!*) 2 - (m av 


±TO rct ) 2 



F F 



i+ v2 )(pi ct - 



(1 



4(i 



(B12) 



A similar expression can be obtained from (|B5|) for lljj 6 ; . 

For T = 0, the integrand in (|BT2|) can be simplified. Using JBTJ and £* et (fc) = E av (fc), we find ImE c (fc) = 
sgn(fco — /i)Im E ret (fc) = — sgn(fco — /i)Im E av (fc). This means that the time-ordered imaginary parts are identical to 
the retarded ones for fco > [i and identical to the advanced ones for fco < M- The real parts of $] c . rot > av are identical 
for all fco, cf. (jBip . In (|B12j) every retarded quantity is connected to a factor (1 — np) while every advanced quantity 
comes with a factor np. It follows that all retarded and advanced quantities can be replaced by time-ordered ones. 
After the replacement ret, av — > c, the four terms in the integrand can be added up and we find 



ns,(*) 



T=0 



2iN f N, 



*P 



(P+ — P-f — (™+ i m c _f 



(2tt) 4 (fc 2 -m c2 )(p c _ 2 



(B13) 



Before we continue, it is important to realize that energy independent terms which are not covered by dispersive 
integrals can be isolated in the fco — > ±oo limit. By inserting the effective masses and momenta explicitly, we find for 
the numerators on the rhs. of (|B13[) in the a and the 7r case 



(p+ - p-f - (m+ + m_) 2 = fc 2 -4A/ 2 - 2fc M A E M +A E M A S M , 
(p + -p- f - (m+ - fh-f = fc 2 — A E 2 - 



-2fc M AE M +AE M AEr 



(B14) 



with A E 5,A1 = — E^^ 1 and M = mo + ^(^+ + El). The numerators depend on po and p only in the arguments 
of the self energies. Up to the A E terms, the expressions in (|B14|) resemble the mean field results fc 2 — 4m* 2 and fc 2 , 
cf. (|A7IA4I) . 

At large fc 2 , the imaginary parts as well as the four- momentum dependent, dispersive real parts of the self energy 
will vanish due to the cutoff of the model. Only the mean field self energy and the three- momentum dependent, real 



K S J\P±\). Here E* ff is the 



shifts to the dispersion integrals remain finite, i.e., E^_ — > E mf + Ep ffs (|p±|) and E^_ 
effective Fock self energy that is defined below in Eq. (|C5[) . 

The po dependence of (|B14[) is completely lost for fc 2 . — > oo: The fco dependent part of the numerators reduces 
to fc 2 - 2fc M AE M . The term 2fc M AE /i - that is not present in the mean field approaches - does not contribute to 
the integral in (|B13[) . The p + and the p- contributions in AS cancel each other at large fc§. Hence, the energy 
dependence of the numerators in (|B14[) is - as in the mean field case - given by 



(p + — p-f — (m + ± m_ 



fc 2 



O(k° ). 



(B15) 



The same result can be found for the more general expression in (|B12[) . At large fc 2 . , where the self energies are purely 
real, the distinction between E ret and E av becomes irrelevant in the numerators. Note that the second term on the 
rhs. of (|B15p . O(fco), is still p± dependent. We will come back to that point in Appendix [Cl 

Let us briefly comment on the fc — > behavior of the numerators in (|B12IB13[) . This is important for the discussion 
of the pion mass in Section IIII El The integrand of (|B13|) vanishes for fc — > in the pion case since p c + — p c _ and 
m c + — rh°_ for fc — > 0. In (|B12|) the terms that contain either only retarded or only advanced self energies will 
vanish in the same way. The mixed terms become zero already for fco = 0, due to the distribution functions: When 
Pq = Pq = po, the factors (1 — np )np : could only be finite for po < fi and po > fJ, at the same time. 
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We turn our attention now to the denominator of (|B13|) and introduce the integral 

= / (27^4 (pC + 2 _ ^c + 2 )(?5 c_2 _ ^c_ 2) • (B16) 

Note that the imaginary part of I corresponds to the real part of n c,rct . We can introduce a retarded integral via 
ReI rot (fc) = ReI(fc)[l + 2n B (fco)]7 Iml rct (fc) = Iml(fc), cf. (|B1|) . In the mean field approximations, Rel(fc) and Iml(fc) can 
be calculated analytically when a simple cutoff scheme is used. On that level, Iml ret (fc) calculated from a dispersion 
integral is identical to ImI rot (fc) obtained from the direct calculation, cf. Appendix [A"l No constant shifts must be 
added to the dispersion integral (|A8|) . 

The complex self energies that enter (|B16j) in the 0(1/N C ) approach do not introduce additional poles to the 
integrand. They remain small and vanish or become real and energy independent for energies a few times larger than 
the cutoff A. Thus, the full calculation does not introduce additional shifts to the dispersion relation and - like in the 
mean field case - the relation (| A8|) should hold. The complex integral I ret (fc) of the 0(l/N c ) approach is found by 
inserting the second term of ()B8|) into (|B5[) . cf. Eq. (|B12[) . and replacing the numerators of all terms in the integrand 
by 1, 



I ret (fc) = 



d 4 P 

(2tt)4 



' (%-n+) (B17) 



(p T f 2 - m r f 2 )(p^ 2 - m^ 2 ) 
1 



2\ F f £ew2 ^.av2\^^;av2 _ ~ av2\ 



Using V and W from (|16|17|l we find for the real part of I ret explicitly 

Rel ret (fc) = / -7^~—~rr - — n — 2 ^t W 9~ 9-r( n F - 4) • (B18) 

w J (2tt) 4 (v 2 + + wl)(v 2 _+wly F F 

Eq. (|A8|) does not translate directly to the polarizations. The numerators (|B14[) in the integral for n^ e ;(fc) are more 
complicated than the trivial ones in (|B17[) . A factor k 2 occurs which changes the k 2 — > oo behavior. This factor will 
generate a fcg independent shift when it is included in the dispersion integral, see (|A9I) . 



fc 2 Iml( ret >(fc , k) = - f°° dp 2 ReV ct ( PQ ,k) + -V f°° dp 2 ^ k 2 )Rel^(p Q ,k) 

7T JO * JO k 0~ PO 

The numerator in the last term on the rhs. is just p 2 ReI rct (p) for p = k, i.e., this term is a dispersion for fc 2 I(fc) instead 
of I(fc) alone. On account of (|B14|) this is just one of the terms which appear in (|B13|) . 

The other terms in the numerators of (|B13p should not generate additional shifts: They are either constant or - as 
discussed earlier - well behaved functions that vanish for large kg . Hence, the integrals 



d 4 p 


4Af 2 4 


2fc Al AE M - 




(2tt)< 


{P% 2 


— fn c + 2 )(p c _ 


2 - m c _ 2 ) 


d A p 


AS 2 4 




AE )1 AE' J 


(2tt)< 


(Pi 2 


— m c + 2 ){p c _ 


2 - ml 2 ) 



(B20) 



cf. (|B14j) . will satisfy corresponding dispersion relations without shifts (|A8|1 as I(fc). The integrals i r °\(k) can be 
constructed by inserting the second term of (|B8ll into (|B5|) . bringing the numerators into the form (|B14j) . and 
removing the k 2 terms. 

can be decomposed into 11^ (A;) = 2iN f N c [k 2 ¥ ct (k) - J^„(k)]. The dispersion relation for U^J is then given 
by the dispersion relations for fc 2 I ret (fc) and 3 I ° t 7r (k). A shift proportional to the one in (|B19|) will be the only correction 
to the dispersion integral over Imll™;. Using (f30|) we find 



Rengjft) = lv rdp f^f ~ 2JMl I dp 2 Rel-\ P0 ,k). (B-l. 

ft o Po 



Eventually, the dispersion relation for the retarded RPA polarizations - including all fco independent shifts - is found 
by combining the results from Eqs. (|B11[) and (|B21[) . The real parts which enter our numerical calculations are 

Re II[ ct (Jfe) = n n +Re (k) . (B22) 
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APPENDIX C: REAL PARTS OF THE QUARK SELF ENERGY 

The 0(1/N C ) meson exchange diagram in the first line of Fig. [4] may contain ko independent, non-dispersive 
contributions to the real part of the quark self energy. Searching these components is not as important as for the 
mesons. The quark mass is determined predominantly by the 0(1) mean field self energy £ mf , cf. Section IV Gl 
that we calculate separately. In the following, we will use the relations between the different kinds of propagators 
(S c ' ret ^, A c ' ret '^) from Sections HID] and HH] without further reference. 

Using Feynman rules, we find for the time-ordered quark self energy 

-»E c (fc) = ^S mf + ^ J ^T l S c (p)f l A c l (p-k). (CI) 

An investigation of the integral in the same way as in Appendix [Bl is not possible: In Appendix [Bl the explicit structure 
of the quark propagator is utilized to decompose and simplify the integrand. The structure of the RPA propagator 

Of 

A ^ = - i + 2G n ?(fc ) <™ 

is completely different from the quark propagator. Thus, the methods of Appendix [B] are not applicable here. It is 
not feasible to replace A[ by a standard meson propagator, D^T 1 (fc) = k 2 — m 2 — H c (fc) , in (|C1[) since the properties 
of and differ - especially in the high-energy properties. The propagator D c will vanish in the k 2 — > oo limit. 
In our effective propagators, the polarization should be - up to \k | independent terms - suppressed by the cutoff. 
Therefore, A c will approach a finite value, cf. (|C2j) . 

We split up the RPA propagator into an energy independent part A°°(|fc |) that corresponds to the limit at high 
k 2 and an effective propagator A^ that will vanish for k 2 — > oo, 

A?(fc) = A~(|£|) + £?(*). (C3) 

The same decomposition can be done for A[ et ' av . Since A°° is real, we find £z(fc) = — 2ImA[ et (fc) and, using (|23|) . 
A*(*) =£*(*). 

Inserting the decomposed RPA propagator into the time-ordered self energy (|C1|) yields 



£ c (£0-£ mf = £| -0^ [ir ; s G (p)r ; A°°(^-fc|) + ir ; s c (p)f ; Af(p-fe)] . (C4) 

The first term is a fco independent contribution to Re£ c . Thus, it is not covered by the dispersion integral. In the 
simplest case, A°° = — 2G (for IF — > 0), this term becomes just the constant Fock self energy. This is not surprising 
since the Fock diagram is part of the effective meson exchange, cf. Fig. [3] 

Because of its structure, we will refer to the first term on the rhs. of (|(J4|) as effective Fock self energy E^ ff . Using 
the decomposition S c (p) = ReS ret (p) - |-A(p)[l - 2n F (p )] and Eqs. (fTUITTj) . we find 

E \ s (|fc|)=2| -^A°°(\p-k\)A s (p)MPo), 

^U\k\)=2j ^A~(|jf-k|)^(p)[l-2n F (po)], (C5) 
S e F ff ,,(|fc|)=-4 J ^A°°(|p-fc|)cos^(p)n F (po). 

Note that Ao is symmetric in po in the vacuum while [1 — 2np(po)] is antisymmetric. Thus, E^ ff becomes finite only 

in the presence of a medium. Since A v (p) = A v (po, \p\) and J_ 1 (icos^cos i? = 0, Sg ff „ will be zero when A°° is 
constant. 

As discussed before, a search for shifts to the dispersion integral in the second term on the rhs. of Eq. (IC4I) is 
complicated. We just give a plausibility argument here, why such shifts should not occur. In general, the existence of 
constant shifts is related to the level of divergence of the integrand in the dispersion integral. The RPA polarizations 
(|Cip include two quark propagators ~ P±/p± (for large k 2 ). This gives rise to a term ~ k 2 in the numerator of 
the integrand. We have shown in Eq. (|B19p how the inclusion of the factor k 2 in the dispersion integral for Rellj ct 
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generates a constant term. The other terms in the numerator of (IC1[) are of lower order in k and do not generate such 
shifts. 

In contrast to the polarizations, the self energy integral in (|C4j) contains only one quark propagator and one effective 
meson propagator A 1 . Assuming that A^(p± k) ~ l/(p± k) 2 for large k 2 - like for the standard propagator - the 
last integrand of (|C4[) does not contain a factor k 2 in the numerator. The denominator, on the other hand, resembles 
that of the polarizations. Consequently, the imaginary part of the quark self energy (IC4|) is of lower order in k than 
the imaginary parts of the RPA polarizations. We can expect that the dispersion integral for Re E rct does not include 
a shift similar to the one in (|B19p . If we overlook any (next-to-leading order) shifts by applying this argument, they 
should be small compared to the leading-order Hartree self energy. 

The complete dispersion relation for the real part of the retarded quark self energy is then given by 

ReS rct (fc ,fc)-E mf + S^) + ^ ( ™ dp Q ^^-. (C6) 

2tt J-oo fc - p a 

The Lorentz components of Re S ret can be found in the usual way ((5]). Note that all components depend on \k | only, 
not on the full three-momentum k. 

A remarkable interplay between quarks and mesons can be observed here. For k 2 — > oo, ReS rct approaches the 
|fc| dependent, finite value of S mf + S^ ff (fc). The momentum dependence is introduced by the high-energy limit of 
the RPA propagators, cf. (|C3IC5p . The momentum dependent limit of the meson propagators, on the other hand, is 

generated by £g ff (fc) that enters Eqs. (|B14IB15[) . This self-consistent effect is not present in the Hartree(-Fock)+RPA 
approximations since the polarizations are not fed back into the quark self energy. There, the k 2 , — > oo limits of Re E ret 
and Re II rct are constant in k. 
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